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We analyze a class of linear wave equations for odd half spin that have a well posed initial value 
problem. We demonstrate consistency of the equations in curved space-times. They generalize the 
Weyl neutrino equation. We show that there exists an cissociated invariant exact set of spinor fields 
indicating that the characteristic initial value problem on a null cone is formally solvable, even 
for the system coupled to general relativity. We derive the general analytic solution in fiat space 
by means of Fourier transforms. Finally, we present a twistor contour integral description for the 
general analytic solution and assemble a representation of the group 0(4, 4) on the solution space. 



I. INTRODUCTION 



It is a well known fact that many spinor equations that are perfectly well behaved in (flat) Minkowski space can not 
be translated to a general four-dimensional curved background manifold. This happens e.g., for the zero rest-mass 
zrm) fields with spin s > 1 and for the twistor equation. In these cases the appearance of the Buchdahl conditions ( 



12 1) imposes algebraic conditions relating any solution of the field equations to the (conformal) curvature of the 
manifold. In the zrm case, these conditions are very restrictive in that they limit the solution space of the equations 
for a general given background manifold. In the case of the twistor equation one finds that solutions can exist only 
on algebraically special manifolds of type N or O. Recently there has been some interest in the case s = 3/2, 
the Rarita-Schwinger system, for which the consistency condition is just Ricci-flatness of the manifold such that one 
can regard the vacuum Einstein equations as integrability conditions for this system of spinor equations. 

Despite this very interesting approach to the vacuum Einstein equations and its relation to twistor theory we want 
to present here a class of spinor equations which does not have the drawback of being well defined only on flat space- 
time (cf. also [^). In particular, these are linear equations for a spinor field of half integer spin s, which include the 
Weyl neutrino equation as the case s = (1/2). The general equation we shall consider may be written as follows: 

= 0. (1.1) 

Here the operator Daa' is the standard Levi-Civita spin connection and the spinor field <I>^b '^cd totally symmetric 
and has m + 1 unprimed indices and n primed indices, with m and n non-negative integers. A simple count shows 
that the field 4>^b ''^cd has (m -I- 2)(n -I- 1) components at every point, whereas the number of field equations is 
(m -|- l)(n -f 2). So the excess number of field equations vis a vis components is m — n. In the case m < n, we have 
fewer field equations than components. Typically this leads to gauge freedom (at least in flat space): for example in 
the case (m,ri) — (0, 1), equation ( |1.1| ) gives the self-dual Maxwell equations for a potential (j)a which has the gauge 
freedom (l)a '-^ (f'a + dacj) with </> an arbitrary scalar field. Next in the case m > n, we have an over determined system. 
This leads to integrability conditions and possible inconsistencies in a general curved spacetime. The classic example 



is the case n = 0, in which case equation (1.1) becomes just the standard zrm for spin (m-|- l)/2, which is inconsistent 
in general, at least for m > 4. The subject of this work is the case m = n, where there are exactly as many equations 
as there are field components, so one might expect that there are no non-progagating degrees of freedom and no 
constraints. Indeed we will show that these equations possess the following properties: 

- They can be derived from a variational principle; 

- They are conformally invariant; 

- The Cauchy problem is well posed in flat and in an arbitrary curved space; 

- There exists an equivalent exact set of spinor fields JlOt for the fields when propagating on a curved background 
and also when they are coupled to gravity, which means that the characteristic initial value problem is formally 
well posed; 

- The fields propagate along null hypersurfaces in flat space. However, they do not satisfy the wave equation 
□0 = but an equation □™+^(/) = instead, where □ is the d'Alembert operator corresponding to the spin 
connection; 
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- The general solution in Minkowski space can be given using a variation of Fourier transforms; 

- A twistor description for analytic flat space solutions can be given. 

The plan of the paper is as follows: in section 2. we present the variational principle and derive the energy momentum 
tensor. To do this we need to explain and extend a formalism given elsewhere In section 3. we discuss the Cauchy 
problem and show the existence of solutions. In section 4. we explain the relation with exact sets and discuss the 
characteristic initial value problem. Section 5. is devoted to the general solution in Minkowski space and in section 6. 
we show how to obtain solutions by performing contour integrals in twistor space thereby establishing an isomorphism 
between analytic flat space solutions and certain cohomology groups for suitable domains in projective twistor space. 



II. THE VARIATIONAL PRINCIPLE AND THE ENERGY MOMENTUM TENSOR 



In this section we will derive the spinor equations from a variational point of view. We will use a formalism described 
in 1^ to derive our results. In order to introduce the notation and also to transcribe the various formulae to apply to 
spinors we will briefly review its basic features. 

The starting point is an SO {1,3) principal bundle 0{M) over space-time M carrying a tensorial one- form 6'° with 
values in and a connection form 0°-t with values in the Lie algebra so(l,3) of the structure group. There also 
exists a constant matrix rjab of signature (1,3) that is used to construct the Lorentz metric on M. We will require 
that the connection be torsion free, i.e., DO"^ = 0; D is the exterior covariant derivative. Making use of the standard 
2— to— 1 epimorphism of SL(2, C) onto the orthochronous Lorentz group to enlarge the structure group of the bundle 
and employing abstract index notation we may write 9"^ = , 9'^b = —e^B9^B' ~~ (-^ B'9^ b, thus defining the 
unprimed and primed spin connections, symmetric in their respective indices (the sign is chosen in order to conform 
with other references). Then the torsion free condition is (note, that we suppress the wedge because it will be the 
only product we use between forms) 

Q^d9^^' +9^b0''^' +9^' B'O^^' . (2.1) 

We introduce a set of vector fields Oaa', dual to the forms 0^^ and 9^b- Their action will be extended from 
functions to indexed quantities by requiring that they be derivations of the algebra of indexed forms annihilating 9^^ 
and 9^B^ see for further details. 

We define a variation at 9°^ and 9°'i, as the derivative at A = of one-parameter families 9'^{X) and ^?''f,(A) with 
6*^(0) = 6*" and 9"'i,{0) = 9°'i,. Denoting the variations by a-nd x"'b we find that x° is a tensorial one-form with 
values in R**, is a tensorial one-form with values in so(l, 3). Using spinor indices and decomposing into irreducible 
parts we have 

X^^' = a^B-^ + ^<j9^^' + T^Be^""' + f^'s.e^^'. (2.2) 

The variations are not independent; in fact, the torsion free condition entirely fixes the variation b — bcc'O'^^ 
of the connection for given x°: 

XABCC = ~d(^A'^B)CE'C' + J^C{BdA)C''^ + ^cc'Tab- (2.3) 

Spinor fields can now be considered as tensorial functions with values in the appropriate representation of S'L(2,C). 
In particular we consider totally symmetric fields "^s^ "c^ with m primed and m+l unprimed indices. Note however, 
that we will switch back and forth between viewing the fields as tensorial functions on the bundle and as spinor fields 
on space-time M (section in an appropriate associated bundle). With these preparations we can now write down the 
following horizontal four-form 



Im 



[^j;,'^c-^'Di.^?::c^,)^AA', (2.4) 



with Tia = {l/^)^abcdO O'^S ■ Notc that we consider our spinors to be commuting quantities. If instead we w ould need 



to have anticommuting spinors we would use the real part instead of the imaginary part in equation (2.4). Since C 
can be considered as the puUback of a unique globally defined four-form on M. we may define the action A as the 
integral of C over Ai: A = /^'C. It is then easily verified that the variations of A with respect to ^/i (■(/') give the 
equation (and its complex conjugate): 
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Dibfj^,:::2"^A')A = o. (2.5) 

Now D = 0"-da and e^Y.b = S^Y. with E = {l/2A)eabcdd''d^9''e'^. Using this and stripping off the form E, we obtain 
the desired field equation 

5a(A'V'b^co = 0- (2.6) 

Let us now discuss some of the basic properties of this field equation. First, we note that for m — this is just the 
Weyl neutrino equation Oaa''^'^ = 0, which is the zrm field equation for spin (1/2). Just as the neutrino equation the 
general equations are conformally invariant if we assume a transformation of the fields with conformal weight —2. We 
define the other irreducible parts of the covariant derivative of ip, 

'^^■■■^ (2.9) 



(2.10) 



'^B'...C' — '-'E VE'A'...C' 

Then we have the expansion 

C)E ,AB...C _ ^EAB...C , ,,EAB...C , E{A^B...C) C9 1 1 ^ 

Oe'Vb'...C' — ^E'B'...C' —7^E'{B'I^...C') —^^E'{B'f^ ^...C) ' 

By virtue of the field equation these fields satisfy the following relations (among others): 

a ,,AC...D n j,ABC...D 

m+ i "MA'^C'...D') =^abVa'C'...d'^ (^-1^) 

™ f}iA BC...D) _ q(A i BC...D)E _ 1 □ , ASC.U 1 3) 

^j^l"(A'^C'...D') EVa'C'...D' 2^Va'C'...D'i l^'-L'^j 

dA(A'\i^.::i^>^>) = □(A'B'V'gFvz^'), (2.14) 

c^A' XAB...CD _ 1 rn ,ac...g p (B,C...I?)A 

'^A'C'.-.D' — ~2(m+ i)^c''---D' ^ VC'...D' + l^-J^oj 

Q ,,AB...CD 1 |-|„/.-B---g-P I n-D' ,/,-BC...U , /'o 1 7^ 

(JA{A'fJ'B'...C") = 2(m + 2) '^('4'V^i3'...C')D' + (2.17) 

™+l|-^ ,C...D)A , 1 a(B C...D) /oio\ 

-^□^(^^^,^,^.^, + j,^, (2.18) 

a| M^?c^.:.g^ = {a^'^^-^;c' ^.z.' + n^'^Vi'^c^^z,'} ■ (2.19) 

In these formulae we have used the spinor curvature derivations Dab and Da'b' as defined in |]l2| . To further analyze 
the situation it is very convenient to introduce four differential operators L, M, M' and N acting on irreducible spinor 
fields by taking the covariant derivative and then projecting onto one of the four possible irreducible parts (see |^ 
for a rigorous definition and further details). Thus, for a field cf) with p unprimed and p' primed indices and all 
its indices down, we identify L(j) = d(E'{E^A...B)C' ...D'), M'(f> = -p'd^^^pA...B)C" ...E' , = ~pdf^,^c' ...D')a...e, 
N(j) = pp'd^^ '4'A...EC' ...E' ■ These operators obey certain commutation rules, most of which are trivial in flat space. 
The nontrivial ones are [L, iV] (/> ~ —\{p+p' + 2)00 and [M, M'\ (j) = — i(p — p')n0. There is a further relation: 
LN(f) — M M'(j) — \p{p' + 1)D0 = 0. The wave operator commutes with all the derivative operators. In flat space, the 
equations ( |2.14| ), ( |2.12| ), (|t|) and ( ^.13D above can be written as follows: 



M'A = 0, MV = 0, M'ii=\UiP, N\ = -m{m+\)Lv+U2m + ?>)U%l^, Liy = -^^^^-^D^;. (2.20) 

^ ^ 2™ 

We observe that v and A are spinor fields of the same class as ip obeying the same equation. In contrast to the zrm 
case, the field V' does not obey the wave equation □■0 — (unless to = 0, because then /i = 0). However, with these 
preparations it is now easy to prove the 
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Proposit ion II. 1 Given a smooth spinor field ip with m primed and m + 1 unprimed indices subject to the field 
equation then O'^'^^tp = in flat space. 

Proof: We prove this by induction on m. The case m = is the Weyl neutrino equation for which the assertion is 
true. Now assume it is true for (m - 1), then = 0. But then = LD"'v = W'^LN^j = □"(^(to + l)^^)^^ = 

i(m+ l)2n"+V. 

FinaUy, we want to derive the energy momentum tensor of these spinor fields. This is usually done by considering 
the action as depending on the metric of the background manifold and then varying with respect to that metric. 
The result is the natural object that would appear on the right hand side of the Einstein equations if the system 
were coupled to gravity. In our case we can not regard the action as depending on the metric, we have to take it 
depending on the canonical one-form. Then the variation of the action with respect to 0" contains terms proportional 
to OABA'B', TAB and cr. The functional derivative of the action with respect to oaba'B' is the trace free part of the 
energy momentum tensor while taking the functional derivative with respect to a gives the trace part. However, in 
the present case, we expect this term to vanish due to the conformal invariance of the equation and the functional 
derivative with respect to tab will be seen to vanish as well. This is related to the fact that the connection is required 
to be torsion free. The variation of the action with respect to 9^^ is 5A — Im J 5C with 

5C ^ {(m + 1)^^:^-^ X^^i.V^frS''' - } ^AA' + {^^ ^-^'i^^^^c^} S^AA' (2.21) 

Using a formula from Q| (which however contains a misprint) we find SJla = 2x''ib'^a]- If we now put all the cr-terms 
in B equal to zero retaining only the t terms we get 

6C^{{m+ l)^^:^^-'^'i?r(^^^f - m{p^[^^-^'Df^\s'i^^P::^g + :.^-^'i?^|/^;.#r-4-^'^^, } E^^,. (2.22) 

Integrating by parts and using the field equation several times gives 

SC = ~i,i[''c-^'T^EDi,iP::cf^AA' + m4,j;,^'c-'''f^' B'D^^'l^v^^AA' (2.23) 

+ i,^^:c--^'f'''A'Di^i?-J,i:AE' (2.24) 

Now the first and third term cancel while the second and fourth term combine to a multiple of the field equation. 
Hence the functional derivative of the action with respect to the t terms vanishes. By a similar argument one can 
show that the terms proportional to a also vanish so that one has to consider only the trace free parts proportional 
to uaba'B'- In this case the calculation is similar but more complicated, so we only state the result. The energy 
momentum tensor of the spinor fields subject to equation ( p^.l| ) is 



T^"-^'"' - Im{(2rn- l)^g:.:,"'^^'Ag:)^^,^-" - "^(^"^+ ' ■■"'a^^^.^v^ 

"^'''"'"'\^^:T-"'^"4^."."j. (2.25) 



(m + l)(m + 2) 

By construction, it is divergence free and due to the conformal invariance it is also trace free. Note, that it is made 
up of the fields and all the non vanishing irreducible parts of its first derivative. The case m = agrees with the 
conventional energy momentum tensor for the Weyl equation. 



III. THE CAUCHY PROBLEM 



In this section we will prove that equation (1.1) has a well-posed Cauchy problem, i.e., we will show that given 
appropriate Cauchy data on a spatial hypersurface S there will exist a unique solution of equation (1.1) on a small 
enough neighbourhood of S. So existence and uniqueness will hol d (o nly) locally in time. 



We will first examine the hyperbolicity properties of equation (1.1). Let us write the field equation in the form 



\d' 



• °B°c) 



We abbreviate the product of (5's by thus introducing the clumped indices a 

...C 

.c 



(3.1) 

AA' and p,, indicating elements 



of the spin space S^^' qI and its complex conjugate dual space. Then equation (3.1) has the form 



(3.2) 



For each covector pa, PaA'^^ defines a map P from S*" into S'p which is easily seen to define a sesquilinear form on 
S'^. Before proceeding further, we will prove two useful lemmas concerning the map P. 
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Lemma III.l The map P is an anti-isomorphism if and only if Pa is not a null vector. 



Proof: For the sufficiency we use induction on the number m of primed indices. For m — the map is i— *■ pa'AI^^- 
If pA'Ai^^ = then p^pA'Ai^^ = —\p^vd = 0, where p^ = PaP"^- Hence, if 7^ then P is injective and therefore 
bijective. Now suppose the statement is true for an integer m — 1; we will show that this implies that it is also true 
for TO. In this case the map is 1— > Pa{A' ^b^"c')- Pa{A'^b^"c') — ^ then also = bPa{A' ^b^"c') — 
:^;^Pa(b'^'^C')A'P^' B- The induction hypothesis implies that p^i-^'A^? -^ = and therefore A' A^^' ^; = but this 
implies p^ ^Paa'^^? "c' = \p^^^?"c' — ^- '^^ 7^ map is injective and therefore bijective. If = then 
we may write paa' = PAPA' for some spinor pA and its complex conjugate and then Xb^' ^' = P^P^ ■ ■ ■p'^'pb' ■ ■ -Pc 
is a nonvanishing spinor that is mapped to zero. 

N 

Lemma III. 2 The determinant of the sesquilinear form defined by P is det(paAp^) = c{paP'^) ^ where N = (to + 
1)(to + 2) is the dimension of S'^ and c is a non-zero real number depending on the choice of basis. 

Proof: Consider the "characteristic polynomial" Q{p) = det(poAp^). As a determinant it is a Lorentz scalar and since 
the only available scalar for a given pa is PaP"' it follows that Q{p) is a function of PaP"'- Since Q{p) is a homogeneous 
polynomial in pa of degree N the result follows. Another way (which is useful later) to see this is the following: The 
determinant is proportional to e'^^ ' '^"Pa^^j,yj • • ■Pa^p„,y„e''^ ' '^" where g'^i - '^w jg (dual) volume form on S'^. As 
such it is built up from the volume form eAs of spin space. After contracting away all the e's in the expression we 
are left with A'' pa's which are all contracted with each other. Due to the identity paa'P^^ = ^PaP^^A'^ each pair 
of Pa's contributes one factor Pap" to the determina nt. Since there are N/2 pairs, the result follows. 



From these two lemmas, we see that the system (3^) is not symmetric hyperbolic unless m = 0. If this were the 
case, then there would exist a timelike future pointing covector pa such that paA°^^ was hermitian and positive definite. 
Although the system is symmetric, it is not definite. For let Pa be any timelike future pointing covector with Pap° = 2 
and choose a spin frame {o^,l^) such that pa = oaoa' + la^A'- Then we choose \^ \c' = 0^0^ . . .o^ob' ■ ■ -oc 
and find that A'^Pa^^j^A'' = if to > 0. This, of course, implies that not all (real and nonvanishing) eigenvalues of 
PaAp^ can be of strictly one sign. In the Weyl case it is well known that the equation can be written in a symmetric 
hyperbolic way. 



To proceed further we determine the characteristics of the system (3.2). These are surfaces locally described by 
the vanishing of a function cj) such that it is not possible to determine the outward derivatives of a function from 
given Cauchy data on th e surface. Hence, on these surfaces the "characteristic equation" det(A? 9q(/)) = holds. 



Because of lemma ( III.2| ) each characteristic surface is a null surface. At each point of M the normal characteristic 



cone defined by Q{p) — coincides with the null cone at that point. However, unless to = the characteristic cone 
has multiple sheets which implies that the system is not strictly hyperbolic. The theory for non-strictly hyperbolic 
differential operators is not as well developped as for strictly hyperbolic or symmetric hyperbolic operators. However, 
in our case we can apply a theorem of Leray and Ohya [Q on non-strictly hyperbolic systems of partial differential 
equations. Their main assumption is that the characteristic determinant Q{p) factorizes such that each factor is a 
strictly hyperbolic polynomia]|^ which certainly is the case here because PaP" is a strictly hyperbolic polynomial of 
degree 2. They show that for Cauchy data on an initial nowhere characteristic surface S which belong to a Gevrey 
class of functions with index a the system has a unique solution in that class. This solution admits a domain of 
dependence, i.e., the value of the solution at a point depends only on the Cauchy data in the past of that point. 

Let S* be a spacelike hypersurface in M and t"' a timelike vector field on M. Define a time function t on M by 
the requirement that t = on S and that f^dat — 1. We define spacelike surfaces St as the surfaces of constant t. 
Given coordinates (x^,a;^,x'^) on S we can continue them off S along i° by Lie transport, i.e., by requiring that the 
lines X* — const, are the integral curves of t"'. Thus, we obtain a coordinate system (t — x'^,x^,x'^,x^) on an open 
submanifold S of M that is topologically S xR. We choose a spin frame (o'^, l-^) such that o"^o"^ — \/2n°, the 

unit normal to St. Then 4'ab...cb' ...c ha-s components ip" = 4'kk' with a = (to + l)fc + fc', < fc < to + 1, < fc' < to 
(such that < a < (to + 1)(to + 2)) and 

/ / A B C D A' B' C" D' fo o\ 

Wkk' ^WA...BC...DA'...B'C'...D' Q_.^OL ...L O . . . O > ...i . (3.3) 



'^A polynomial Q{^) of degree n is strictly hyperbolic iff the cone C = : Q(^) ~ 0} has a nonempty interior such that each 
line through an interior point not including ^ = intersects C in exactly n distinct real points. 
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By taking components of equation ( |l.]| ) we obtain a system of equations of the form {Da denoting the partial derivative 
with respect to x"") 

G%''[t,x)Da^'' + V''p[t,x)i;'' = 0. (3.4) 

Here the functions G^"(i,a;) are functions on M which are algebraic expressions in the metric components gab with 
respect to the coordinates {t,x) — x^, x^, x^). The functions r^(i,x) are algebraic expressions in the coefficients 
of the spin connection. We define a^(x, D) as the linear differential operator in (3^). Then we consider the following 
Cauchy problem: 

a1{x,D)il;^{x) = 0, 

tl'^\s is given. (3.5) 



It follows from lemma ([II. 2) above that up to operators of lower order than N = [m + l)(m + 2) we 
have det(a^(a:, D)) = □"S" where we define the determinant of noncommuting quantities by the usual formula 
det(a^(a:, D)) = X^ttgSjv ^^S^'-^) "^^(1) • ■ ■ °^Tr(Af)- '-^ wave operator with respect to the metric gab expressed 

in the coordinates {x""). This operator is strictly hyperbolic with respect to the hypersurfaces St- We are now in a 
position to prove the 

Theorem III. 3 Let a he a real number with 1 < a < j^_2 ■ If the metric coefficients are in the Gevrey class 

and if the initial data are in the Gevrey class 72"'' ("S) then in a sufficiently narrow strip T,' = {{t,x) : < t < T} 

around S the Gauchy problem ^3.^ has a unique solution £ ^ '^"''(S), whose support is contained in the domain 
of influence of the support of the initial data. S' = E if ^ ^ 01 < j^r^ ■ 

Proof: This is a straightforward application of theorems of existence and uniqueness in §6 of |^]. We only need 
to determine the various integers needed in the theorem. We associate the integers — 1 with each unknown 
function ■0'^ and the integers n" = with each of the equations such that order(a^(a:, I?)) < m'' — = 1. Then 
m — "Y^pimf^ — n^) = N \s the order of det(a^(x, £))). Each of the factors in the principal part in det(a^(a;, D)) is 
equal to the wave operator such that aj — □, ruj = 2 for (j = 1, . . . , TV) and the number of factors is p = N/2. Since 
in our case r = 1 we need to add to each of the integers m'^ and the same integer N/2 in order to satisfy the chain 
of inequalities 

0<r<p<n<n°'<fi,n< ,p < m, (3-6) 

as required in With our choice of the integers we have to'3 = 1 + {N/2), n" ^ 0, p = {N/2) — n — n and 
the inequalities are satisfied. According to the theorem the index a of the appropriate Gevrey classes lies in the 
interval 1 < a < j^z:2- coefficients of are assumed to be in 7^^^~"'^'^"-' (E), those of aj are in 7^^^^'^"-*(S) 
for Q < i < N/2 — 1. Since all the factors are the same, aj = □, this implies that the coefficients are in fact in the 

smallest possible class which is 7;^ ^'^"''(E). The coefficients of the operator a have to be in the class j^^^'^°'\t,). 
Taking all this together and remembering how the coefficients in the operators are constructed from the metric this 
implies that the metric coefficients have to be in j^^'^'^"\'Et). Then the conclusion of the theorem implies that the 
solution of the Cauchy problem is in the Gevrey class 7^'^''^'^"^(E). 

According to this theorem there exists a strong correlation between the spin of the field and the degree of smoothness 
of the space-time that admits a solution of the equation. The higher the spin, the "smoother" the space-time has 
to be. The smoothness is controlled by the number of components of the field, N , which depends quadratically on 
the spin m -I- i of the field. We can improve on this relationship somewhat by using a simplification due to Bruhat 
which is based on the observation that if all the minors in det(a^) have a common factor then this factor can be 
ignored which results in a reduction of the number p of factors of det(a^(a;, D)) and therefore in the Gevrey index a. 
To be more precise, we need to prove the 

Lemma III. 4 In the adjoint matrix of paA°^^j all the entries have the factor (p^p°)'"(™+i)/2 in common. 



Proof: As in the proof of lemma (III. 2) this result can be obtained by "index counting". The adjoint matrix is given by 
the expression e.^^'^'"^" Pa^^^ui ■ ■ • Pa N^pNfiv^'"'^ homogeneous of degree N —1 inpa- If we contract over all the 
indices contained in all the e's in the volume forms we are left with an expression that contains A^— 1 = (m-|-l)(m-|-2)— 1 
Pa's, each with one unprimed and one primed spinor index and has (2m + 1) free indices of either kind. Therefore, 
(A^ — 1) — (2m -t- 1) of the PaS are contracted together, resulting in a factor (p^p'^)™('"+i)/2 in each component. 
This lemma allows us to prove the 
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Corollary III. 5 In the statement of the theorem we can extend the range of the Gevrey index atol<a<l + -^. 
More precisely, if the metric coefficients are in the Gevrey class ■j^'^^''^°'\Ti) and if the Cauchy data are in J2°'\S) 
then the Cauchy problem (3.5) has a unique solution ip^ G 7™~''^'^"''(I]') in a sufficiently narrow strip S' around S. 



Proof: We observe that in proving the existence theorem Leray and Ohya use a theo rem for systems with diagonal 
principal part (see §5 of [Q). To apply this theorem one multiplies the system (1.1) with the differential operator 



corresponding to the adjoint matrix of paA? . This renders the principal part of the resulting system diagonal. Due 



to lemma (111.4) it is enough to multiply with the operator (of lower order) obtained from the adjoint matrix by 
dropping the common factor. Then one obtains the result in a straightforward manner by applying the theorem for 
diagonal systems. 

To end this section we want to make several remarks. 

— The case m — (the Weyl neutrino equation) is the strictly hyperbolic case where we can choose a = oo. This 
implies that it is possible to prescribe initial data with only a finite number of continuous derivatives which is 
a known result for strictly hyperbolic systems. 

— It is also worth to mention that in all the cases there exists a domain of influence, a fact which is taken to 
indicate the hyperbolic character of partial differential equations by many authors. 

— The fact that the smoothness of space-time is strongly linked with the spin of the field is an interesting feature 
of this class of equations that is not present in other spinor equations. One has to say, though, that it is not 
known whether this is a necessary consequence since the theorems only provide sufficient conditions for existence 
and uniqueness. 

— The diagonal system used in the proof of the corollary corresponds to the equation = that we 
derived in the flat case in section 2. This equation is distinguished by the fact that it is a linear equation for 
'i'B^ 'c' such that the coefficients are functions not of the connection but of the curvature and its derivatives 
only. 

— The inhomogeneous equation dA{A'4'B^"c') ~ Xa'b^ c '^^^ treated in a straightforward way and one obtains 
existence and uniqueness of solutions in the same Gevrey class as for the homogeneous case provided that the 
right hand side is in an appropriate Gevrey class, see [Q. 



IV. THE FORMAL CHARACTERISTIC INITIAL VALUE PROBLEM 



In this section we want to discuss the formal aspects of the characteristic initial value problem on a null cone for 
this class of spinor equations. Due to the inherent singularity at the vertex of a null cone this problem is very difficult 
to analyze and, in fact, there are no existence results for many partial differential equations appearing in physics, 
most notably the vacuum Einstein equations. So one has to resort to formal methods to obtain at least results about 
the fcasability of existence theorems. A very useful method to achieve this which is adapted to four dimensions is the 
method of exact sets of spinor fields developed by Penrose . It is based on the observation that in Taylor expansions 
of spinor fields around a point it is exactly the totally symmetric derivatives of the field that determine the restriction 
of the field to the null cone of that point (the null datum). Roughly speaking, if a system of field equations for a 
collection of spinor fields has the properties that the totally symmetric derivatives are algebraically independent and 
if they determine algebraically all possible derivatives of the fields then the collection of fields is said to be exact (see 
pO{ , Q for the rigorous definition). 

It has been useful to employ an algebraic formalism based on the four derivative operators L, M, M' , N (already 
mentioned in section 2) which correspond to taking the four possible irreducible components of the covariant derivative 
of an irreducible spinor field. We will not describe the full formalism here because it would take up too much space. 
Instead we will only give a brief summary and refer for further details to j^. The totally symmetric derivatives of a 
spinor field correspond to applying powers of L to the field. We will call an irreducible spinor to be of type (fc, k') if 
and only if it has k unprimed and k' primed indices (irrespective of their position) . Acting on a spinor field of type 
{k, k') the operators L, M, M', N produce fields of respective type {k + 1, fc' + 1), (fc + 1, fc' — 1), (fc — 1, fc' + 1), 
{k — l,k' — 1). We define the operators H and H' by H(j) = kef) and H'(j> — k'cf) for a type {k, fc') field (j>. 

As we have already mentioned in section 2, the commutator of two covariant derivatives induces commutation 
relations between the derivative operators which in general involve the curvature of the manifold and in addition the 
wave operator. The curvature is characterized by three spinor fields ^, <& and A of respective types (4, 0), (2, 2) and 
(0,0). Before we present these relations we need to define an algebraic operation between two spinor fields cj) and x 
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of respective types {p,p') and {q,q')- The only possible way to combine two spinor fields within the class of totally 
symmetric fields in a bilinear way is by contracting over some of the indices and then symmetrizing over the remaining 
lot. This operation is entirely characterized by the numbers of contracted primed and unprimed indices. So we define 
the bilinear pairings Ckk' by the correspondence 

Ckk'{<P,x) - <I>(c,,,d'Ic'''...D'''' XAi...AkB[,...B'^,E'...F')E...F)- (4-1) 

Then the commutation relations between the derivative operators can be given explicitly as 

[L, N] = -{H + 1)T' - {H' + l)T-^iH + H' + 2)0, 
[M, M'] = -{H + 1)T' + [H' + l)T-\{H - H')U, (4.2) 
[L, M] = -{H' + l)S, [L, M'] = -{H + l)S', 
[N, M] ={H + l)U', [N, M'] = [H' + l)U. 

The operators S, T and U and their primed versions are curvature derivations and act on a field (j) of type {p,p') 
according to 

5^ = pCio'(*,0)+p'Coi'(1>», (4.3) 
Tcj, = Pip - 1)^20' + Pp'Ci,, ($, 0) - p{p + 2)Coo' (A, (/.), (4.4) 
t/0 = p{p -l)ip~ 2)C30' (*, 0) + Pip - l)p'C2i' ($, 0). (4.5) 

The action of the primed operators can be inferred from these by formal complex conjugation. There exists an 
additional relation between the derivative operators, the wave operators and the curvature: 

LN ~ MM' = -iH' + 1)T + \HiH' + 1)0. (4.6) 

The formulae describing the action of a derivative operator on a bilinear pairing are quite lengthy and it is not 
necessary for what follows to present them in detail (see ||]). Symbolically they are given by 

OC(0, x) = E g) + a^Cif, Og), (4.7) 

o 

where O is any of the derivative operators L, M', M or N and ai and a2 are rational numbers determined by the 
bilinear product and the type of 4> and %. 

In 1^ we showed that a collection of fields {(j)j} is exact if and only if two conditions are satisfied: 

(i) all the "powers" L^4'j are algebraically independent, 

in) all the "derivatives" of the fields, i.e., all the expressions scfij where s is an arbitrary string of derivative operators 
are algebraically determined by the powers. 

By algebraic independence we mean that there are no relations between the powers involving only the bilinear pairings 
(and possibly the curvature). In the same spirit we mean that the derivatives are determined by exactly such relations 
in terms of the powers. Several examples have been treated in 1^, jl^]. In a certain sense, the powers form a 
complete and independent set of functions that generate the solution space of the equations considered. 

Before we consider the general case we want to study the flat case to find the exact set structure underlying equation 
(1.1). So let V' be a field of type (m + l,m) satisfying the equation M'lj} = 0. Referring again to we see that 



in that case the expressions L'-M'^M'^ N^'^D'-ip for positive integers I, k, j, n and i generate the solution space. By 
use of the field equation and the commutation relations we have j = since all terms with j > van ish (note 
that in flat space □ commutes with every derivative operator). Similarly, i < m because of prop ositi on (II.l) and 



n + k < m because M and N each contract over one primed index. Due to the additional relation (4^) we have the 
following: L^M^mn-i}} - L'M'^nAf"^*- V - L'+iM'^iV^+iQ*- V ~ • ■ • ~ L^+'M^m+^iP where a ~ means "a is 
expressible in terms of 6" . So we find that in the flat case the solution space is generated by the powers of the functions 
iiki = M^'^N^ip, with < fc < m and < i < fc. Note that there are ^(to + l)(m + 2) of those functions. This 
number is in agreement with the general observation that the number of null data per point for a partial differential 
equation is half the number of Cauchy data per point. 

We now claim that the same set of fields is also a generating set in the nonflat case. Before we prove this statement 
we need some more preparation. Let Sn denote any string of derivative operators of length n. We say that Sn is in 
normal order if and only if it has the form s„ = V-M^N^M'^ with l + k + i+ j = n. With each string s„ we can 
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associate a unique normally ordered string of the same length in the following way: if s„ does not contain M' 
then is the unique normally ordered string that contains the same number of operators as s„ does. If there are M' 
operators in s„ we first replace each pair (M, M') which need not be adjacent with (L, N) until there is no M or M' 
left. Then we bring the result into normal order to obtain s^. Thus we get normally ordered strings containing either 
M or M' but not both. Upon applying these normally ordered strings to i/? we get zero for all strings containing M' 
and for those strings with k + i > m. The others result in L^M'^N^ip — L''ipk+i,i', these functions and their complex 
conjugates will be called a normally ordered derivative or a "power". For each power V"4>ki we call Z + fc its order. 
Furthermore, we need to formalize the structure of the relevant terms that will be encountered. 

Definition IV. 1 A t-term ("t" for "tree") is recursively defined either 

(i) as a power L'^ipki or V"ijJki for nonnegative integers I, k, i or 

{ii) to he of the form Ckk'iR,t) where R is any derivative of any of the curvature spinors and their complex conjugates 
and t is a t-term. 

A t-expression is the formal finite sum "Y^- ajtj with coefficients aj € Q and t-terms tj. A t-term that is not a power 
will be called a pure t-term and a pure t-expression is a Q-linear combination of pure t-terms. 

The pure t-terms are binary trees with Ckk' as nodes and with powers or (derivatives of) curvature spinors as leaves. 
In fact, exactly one leaf is a power all other ones are curvature derivatives. This reflects the linearity of the system. 
We are now ready to prove the 

Lemma IV. 2 Let -0 be of type (m -I- l,m) and satisfying the equation M'ip = 0. Let s„ be an arbitrary string of 
length n of derivative operators. Then Snip can be written in a unique way as Sn^jj — as^'tp -\- t where a € Q and t is 
a pure t-expression which contains only powers of order strictly less than n — 1. 

Proof: We use induction on the length of the string. With n = 1 there are four possibilities: Lip — LipQo, Mip — ipi^, 
M'lj} = and Nip = ipii\ so the statement is true. 

Let O denote any of the derivative operators and assume the statement to be true for all strings s of length less 
than or equal to n. Then consider [Osn)ip = 0{snip). By the induction hypothesis and linearity of O w e need to 



consider only two cases, namely Snip is (i) a pure t-term or (ii) a power. In case (i) we need to employ (4.7) to apply 
O to a bilinear pairing, thus bringing O inside the Ckk' to act on each of its arguments. Note, that then O is not 
necessarily the same operator we started with. When it hits the left argument, O converts a curvature derivative into 
a higher one thus producing a t-term of the required type. The other argument is again cither a power or a pure 
t-term. In the latter case we continue descending down the tree structure until we finally hit the power. Then we 
need to consider OL^ipki- By the induction hypothesis this is a derivative of ip of order l-\-k-\-\ < n — 1 and therefore 
equal to the sum of the corresponding normally ordered derivative of ip and a pure t-expression with powers of order 
less than n — 3. The normally ordered derivative is (when non vanishing) equal to a power of order n — 1, hence the 
application of O to a pure t-expression yields a pure t-expression of the required type. 

In case {ii) we have Snip = L''ipki with I -\-k = n. Let us first suppose that I > 1. Then OL'-ipki = [O, L] L'-^^ipki + 
LOL^~^ipki. In the second term we can replace OL^~^'tpki with the sum of the normally ordered derivative and a 
t-expression by the induction hypothesis. Then applying L yields a normally ordered derivative of order n + I and, 
as was just shown, a t-expression of the required type. So we are left with the commutator term. If O = L we are 



done, li O = M ov O = M' the commutator term is equal to a linear combination of curvature terms by (4.2) and 



(4.3)-(4.5) which are t-terms of the required type. When O — N we obt ain apart from curvature terms as before a 



term involving the wave operator. This term can be rewritten using (4.6) as a linear combination of curvature terms 
and the terms LN L^~^ipki and MM' L^~^ipki. The first term has been shown above to be of the correct type and with 
a similar argument one shows that the second term is also. 

Now suppose I = and k > i. Then we need to look at a term of the form OMM'^^^^^N^ip. The only nontrivial 
cases are O = M' and O = N. In the latter case we find that the commutator term is a curvature term and therefore of 
the correct type. The other term is shown to be correct by similar arguments as above using the induction hypothesis. 
In the case O = M' only the wave operator term appearing in the commutator term needs a different treatment. But 
this has been shown above also to lead to correct terms. 

The last case is / = and k = i. Then we are looking at ONN^^^ip. Here all cases are trivial except for O = N 
and this case is treated as above. So, in summary, we have shown that all the appearing terms are of the stated type 
and hence the lemma is proved. 

From this result we can obtain a set of equations satisfied by the functions ipki- Consider M'ipki, a derivative of ip 
of order -I- 1; therefore, there exist equations 

M'1Pk^ = ak^L1Pk^+l + 4 (4.8) 

and similarly 
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M'^|}k^ = V'fc+l,^, (4.9) 
N^l)k^=^Jk+l^^+l+tk^, (4.10) 

where aki = if fc = i and where tki and t'j,^ are pure t-expressions which contain "0;^ with I < k — 1 and possibly 
Lil^ij with I < k ~ 2. If we regard these equations as the field equations for the fields tpki then we can state the 

Theorem IV. 3 A formal solution of ( jj. j|J gives rise to a formal solution of (4^ )-{ J^.li ) and vice versa. The set of 
spinor fields {ipki : < k < m,0 < i < k} is exact. 

Proof: The equivalence of the two systems is obvious. We need to show the exactness. Here, condition (ii) con- 
cerning the completeness of the powers is an immediate consequence of the lemma. The condition (i) concerning the 
independence of the powers can be verified as foll ows. Any relation bet ween the powers has to be generated by the 



application of the commutation relations and (^^) to the field equation (LI) and all its derivatives. From looking at 
the structure of these relations one finds that they can not link any derivatives that contain more than two adjacent 
L's. So all the relations that can be generated must already be conditions on the derivatives of the field equation. 
But this is a condition only on derivatives of the form s„M'-0 and not on any power. The other possible source for 
conditions on powers come from the (derivatives of the) defining equations of the ipki- ipki — M''~'^N''ip. However, 
these are not algebraic relations but differential relations between the functions and — upon taking derivatives — 
between the powers. So there can not be any relations between the powers, which therefore are independent. 



This theorem shows that the characteristic initial value problem for the equation (1.1) is formally well posed. This 
is, of course, a rather weak statement, implying only that one can prescribe certain components of derivatives of tp 
on the null cone of a point in an arbitrary way and that this is just enough information for a unique solution to exist 
on the level of formal power series. 

The exact set {ipki} is not invariant (cf. |jl^) because in the expressions for the derivatives in terms of the powers 
there appear the curvature spinors together with their derivatives which are taken to be known background quantitites. 
Thus, these expressions depend on the actual point in space-time that is the vertex of the null cone. Since the field 
equation comes from a variational principle and since, therefore, there exists an energy momentum te nsor we can 



couple the system via Einstein's equation to the curvature. Thus, we write Gab = 87rTab with Tab from ( 2.25 ). Then 
we know how to express the curvature spinors $ and A in terms of ip and its first derivatives. In fact, A = due to 
the conformal invariance of the equation. We can interpret ip as describing some kind of matter field whose energy 
content creates the curvature of the manifold. We have one more unknown function to consider, the Weyl spinor 
which is subject to the equation (a part of the Bianchi identity) M'^f = 2A/<I>. Referring to a theorem in |^ we 
see that the enlarged set {V'fci,^'} will be an invariant exact set on M provided that we can show that M$ is a 
t-expression and that = 0. We need to interpret the term "t-expression" a little different now because whenever 
A or $ appear in the expressions we need to substitute their resp. representations in terms of the fields ipki. Thus, 
we obtain an actual tree structure built from the bilinear pairings whose leaves consist only of powers L^ipki and L'^* 
and their complex conjugates. This reflects the nonlinear nature of the coupling to gravity. The conditions above are 
easily verified, in fact, = is just the condition that the energy momentum tensor be divergence free and since 
$ itself is a t-expression its derivative is also a t-expression as was shown above. So we have effectively proven the 



Theorem IV. 4 The set {ipki,^} subject to the equations ^^)-(^.lL), Einstein's equation and the Bianchi identity 
is an invariant exact set. 

Thus we can make a similar statement as before concerning the system coupled to gravity. The formal characteristic 
initial value problem is well posed. In this case, we do not have a similar result for the Cauchy problem. 



V. THE GENERAL SOLUTION IN MINKOWSKI SPACE-TIME 



Our aim in this section is to present the general solution of the field equation (1.1) in flat space subject to suitable 
initial and boundary conditions. Since each such solution is also a solution of D"^ip — for some positive integer m 
we will first derive the general solution of that equation. Since this does not depend on the existence of spinors and 
on the dimension of space-time we present the result in a slightly generalized form for arbitrary space-time dimension. 
So we are working in M = M^^lT^- Then we will specialize to four dimensions and restrict the kernel of to those 



spinor fields that do satisfy (1.1). Since in flat space we are dealing with a partial differential equation with constant 
coefficients the general solution could be found using methods from the theory of distributions. We will, however, not 
pursue this here but present a different approach which fits better with the applications we have in mind. 

We begin by introducing certain rings of functions associated to the null cone of momentum space. A function f{k) 
defined on the null cone in "fc" -space will be said to be admissible if and only if it obeys the following conditions: 
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(i) f{k) is defined for all null vectors fc", 



(ii) /(0) = 0, 

(iii) the function / is smooth on the complement of the origin, 

(iv) limt^o {t^ f{ik)) = 0, for any real number r and for any non-zero null vector fc"; this limit must be uniform on 
compact subsets of momentum space. 

Condition (iv) controls both the "infrared" and "ultraviolet" behaviour of the function f{k). Condition (iv) is 
needed to guarantee the differentiability and integrability of Fourier transforms involving the function /(fc). 

Denote by K, and K~ the rings of all admissible functions, all admissible functions that vanish identically on 
the past null cone, all admissible functions that vanish on the future null cone, respectively. Note that we have the 
vector space direct sum decomposition: K = ® K~ . 

For any subspacc i? of a ring and any vector variable X denote by R[X] the space of all polynomials in the vector X 
with coefficients in the subspace R. In particular if R is itself a (sub)-ring, then R[X\ is a ring. For any non-negative 
integer m, denote by .Rmi-'^] the subspace of the space R[X\ consisting of polynomials of degree less than m + 1 
and by R^"^^[X] the subspacc of Rm[X] consisting of all polynomials homogeneous of degree m in the variable X. 
In particular, for x a space-time vector-valued variable, we have that every element (j){x, k) of the ring K[x\ has an 
explicit expression as a polynomial in the variable x of the following form: 

Cl>{x,k) = X'^'X'^- . . .x''^^a,a,...aAk)- (5-1) 

r=0 

Here each coefficient tensor, 4>aia2...a^{k), is completely symmetric and is an indexed element of the ring K. Also only 
a finite number of these coefficient tensors is non-zero. Henceforth, each infinite sum we encounter will have only a 
finite number of non-zero terms. 

Denote by da the derivative with respect to the variable x and by □ = gabd'^d^ the wave operator, regarded as an 
endomorphism of the space K[x]. Denote by L[x] the kernel of this endomorphism and define L„i[x] and 

as the intersections of the space L[x] with the spaces i^mla^] and i4r+[a;], respectively. Consider the operator 

kad" as an endomorphism of K[x]. Since the operators □ and kad'^ commute, kad"' restricts to an endomorphism, 
denoted D, L[x]. Note that D is the derivative operator along the generators of the null cone restricted to solutions 
of the wave equation in K[x]. 



A. The operator D 

Proposition V.l The operator D : L[x] — > L[x] is surjective provided the space-time is at least three-dimensional. 

The proof of this first technical result is rather lengthy and proceeds in several steps. We first perform a decomposition 
into space and time to obtain an expression for a general element of L[x]. First pick a unit timelike future pointing 
vector t"' and denote by S the orthogonal complement in space-time of the vector t". We shall use lower case Latin 
indices from the middle of the alphabet to label the (spatial) tensors of S and shall write for the negative of 
the (flat) metric induced on S from the ambient space-time metric. Then the position vector decomposes as 
2^" = (^iC')i where t = x°'ta and one has the relation x°'Xa = t^ — ^'^j, where the tensor 7jfc and its inverse are used 
for index lowering and raising for spatial tensors. Correspondingly, the operators □ and D decompose as □ = df — A, 
where A = jikd^d^, and D = ndt — {n ■ d^), where fc" = {k, k*), k = k'^ta and k- = Kid^. Note that since the vector 
k" is null, we have the relation = k'Kj. 

Given (f>{x,k) € L[x], define (j)o{£,,k) e K[£_] and ^i{^,k) G K[^] to be the restrictions to the subspace S of the 
functions (/)(x, k) and tad°'<j){x, fc), respectively. 

Lemma V.2 The mapping ps ■ L[x] Kl^]"^, (j) t-^ ((/)o,(/)i) is an isomorphism, mapping each solution of the wave 
equation to its initial data on S. 

Proof: Given the pair ((/>o(C, k), 0i($, k)) G ^^[C]^) the function (j){x, k) = Ps^{{(/>o{C, k)) may be given by the 

following explicit formula: 

0((t, 0, k) = cosh(iAi/2) 0o(^, k) + A-1/2 sinh(tAi/2^ 0i(e, k). (5.2) 
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Here the function s co sh(tt) and u"-'^ sinh(M), with u an operator, are to be interpreted as formal power series. Note 
that in equation ( ^.2|) there are no problems with the square root of the Laplacian, since the functions cosh(w) and 
sinh(u) are both even. Also there are no convergence problems, since the functions and 0i are polynomials in 

the variable 

Define the operator A : by A = psDpg^. Then in view of lemma (V.2) we have to show that A is 

surjective. We derive from equation (5^) the explicit formula for the operator A, valid for any pair ((/)o,(/)i) G K[£]'^: 



A(0o, M = (^01 - (k • d^)(l>o, KA(j)o - [k ■ 

So we must now solve the following pair of equations: 

Kf3 — {k ■ d^)a — 7, 
KAa — {k ■ d^)f3 — S. 



(5.3) 



(5.4) 
(5.5) 



In equations (5^) and ( |S.5|) the pair {j,6) is a given element of the space and the desired solution is the pair 

(a, /?) which must be shown to lie in K\£] '^. N ow it is clear from its definition that if [C]^ is cl osed under multiplication 



or division by k, so we may use equation (5^) to eliminate the function f3 from equation (5^). This gives the following 
equation: 



(A - (n . d^f 



(5.6) 



Here we have put n ■ = riid^^, with rii — k ^Ki, a unit vector and a = k '^{kS + (k • d^)j) G K[^]. Note that the 
desired result is false in two space-time dimensions since the left hand side of equation ( |5.6| ) then van ishes identically, 
but the right hand side need not vanish. So we have reduced the problem to solving equation (5.6), given a € K[S] 
such that the solution a must also lie in the space K[£]. 

Proof of Proposition (5.1): We first prove the proposition for the special case with a of the form 



(5.7) 



Here the numbers p, q and r are non-negative integers and the function Vr G K[S^] is homogeneous of degree r in 
the vector var iable and obeys both of the differential equations Aw^ = and (n • d^)vr = 0. It is easy to solve 
equation ( ^.6[ ) in this case explicitly: a solution is just a = [[q + l)(n + 2r -\-2q — 2))~-'^(^^ — (^ • n)^)i7, as is easily 
checked, by differentiation. This solution clearly lies in the space K[£] and is of the form (^ • n)^'r, where r statisfies 
the equation (n • 9j)t — 0. 

The rest of the proof consists in a demonstration that the general case can be reduced to this special case by 
decomposing a into a sum of appropriate terms and then using linearity of the operator A. We first decompose a as 
a sum of terms as follows: 



r=0 



Each coefficient cr,. is required to obey the differential equation [n ■ d^)ar 
for the quantity a^, valid for any non-negative integer r: 



(5.8) 



0. Explicitly one has the following formula 



E 

s=0 



(-1)^ 



(5.9) 



In particular it is clear from equation (5.9) that each function ar belongs to the space K[£]. Note that the operator 
A — (n • d^)^ commutes with the multiplication operator (^ • n). Also if a G A'[^], then we also have (^ • nY a G A'[^] 
for any non-negative integer r. So using the linearity of equation ( ^.6[) and the decomposition of equation ( ^.8|) and 
(5.9) it suffices to prove the solvability of equation ( ^.6| ) with both the functions cr and a lying in the kernel of the 
operator n ■ d^. Denote this kernel (a subspace of the space K[^) by N[^]. 

Next we use a standard fact from tensor theory that any symmetric tensor may be decomposed into tracefree parts. 
In the present language one has, for any r G iV[^], a decomposition: 



y — 



2\r _ 



(5.10) 
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Here we have put = f^^'. In equation (5.10), the coefRcients must obey the differential equation: 

(a — (n • d^)^) Tr = and must lie in the space iV[^]. Indeed for the case that r is a homogeneous function of 

non-negative integral degree m in the vector variable ^, each function Tr may be given explicitly by the following 
formula: 



Z^v 2'-+2''r(A-r + i)r(s + 1) ) J \ k ui 



(5.11) 



with A = (ri — 4 + 2m)/2. It is easily checked by differentiation that the function of equation ( p.ll| ) lies in the 
kernel of both the operators A — (n ■ d^)^ a nd n ■ as required. The proof of compatibility of equations (5.10) and 
(5.11) follows immediately from the lemma (|V.3| ) given below. 

Since every r e K[^] is uniqu ely a sum of its homogeneous parts and each of its homogeneous parts lies also in 
the space K[^], equations ( ^.10 ) and ( ^.11 ) hold also for inhomogeneous functions t S N[^], provided that equation 
(5.11) is rewritten as follows: 



E 



i-iy 



2\r+s 



(5.12) 



where /i = (n - 4 + 2^a9^)/2. 

Combining these res ults and again using the linearity of equation ( [3.6| ), we see that it is sufficient to prove the 
solvability o f equ ation ( [5.6[ ) with the function a being of the form given in (|]^) above. This completes the proof of 
proposition ( 5.73 ). 

Note that since the above proof is compatible with homogeneity, it also shows that the restriction of the map D to 
the subspace Lrn[a;] has range the subspace L,„_i[a::], for every positive integer m and that the restriction of D to the 
subspace of L*^™-'[a;] is surjective onto the subspace L(™~^)[x]. 

Lemma V.3 For j a non-negative integer, define a function gj{z) of the complex variable z as follows: 



9j 



(z — 2r)r(z — r — i) 



r=0 



T{z-r + l)r(r + l)r(j - r + 1) ■ 



(5.13) 



Then the function gj vanishes identically unless j = and the function go is the constant function with value one. 

Proof: First, the case j = is easily checked by inspection. So henceforth assume, for convenience, that j is a fixed 
positive integer. From its definition it is clear that the function gj{z) is a rational function of the variable z and that 
liniz^oo gj{z) = 0. Therefore, the r esult will follow if it is proved that the function gj{z) is periodic. But one has the 
following relations, using equation (5.13): 



9] 



(^) = E(-i)' 



[z — r)T{z — r — j) 



r=0 



r(z - r + i)r(r + i)r(j - r + l) 



E(-i)' 



V{z- 



-J) 



E(-ir 



r(z — r — j) 



r=0 



r(z-r)r(r + l)r(j-r + l) 



EM) 



r(z-r + l)r(r)r(j-r + l) 
r(z - r - j - 1) 



r=0 



-E(-i) 



r=0 



r jz-r- j - l)T{z - r - j - 1) 
r(z-r)r(r + l)r(j-r + l) 



EM 



r=0 



T{z - r)T{r + l)T{j - r) 

Y U - r)r{z - r - j - 1) 
' T{z - r)T{r + l)r{j - r + 1) 



= \^( 1 r '-^ ~ ~ ^^^^^ -r-j -I) , _ . 

' T{z~r)r{r+mj-r + l) 



(5.14) 



So the function gj{z), for j > Q is periodic and therefore vanishes identically as required. 



B. The Fourier transform operator T. 



Consider the (n — l)-form S defined on the null cone in momentum space defined by the formula = 
dk'^^dk"''' ...dfc""-i = fc_^^gaoaia2...a„_i;=;^ Restricted to the null cone one has kak" = kadk°- = 0, so the left hand 
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side of this equation is orthogonal to the null vector The (n — l)-form S factorizes according to the formula 
k^-E = dk°-uj, where the n - 2-form uj is defined by the formula: k^°-^dk°-^ . . . dfc""-il = fcaoe"o"i"2---a„-i^. Both the 
forms S and uj are closed: dE. — dio = 0. 

Given any A(a:;, k) € K[x], we consider its generalized Fourier transform ^(A), which is a space-time field given by 
the following formula, valid at any point x of M: 

J^[\){x) = I e'^"""' X{x,k)E. (5.15) 



The integration in equation ( 5.15 ) is to be carried out over the complete (past and future) null cone, equipped 
with the induced orientation from its embedding in fc-space, which in turn is oriented by the volume form, 
£ ai a.,.. ..a„ dk°'^dk°''^ . . . dk""" . Note that by definition of the space K[x] the convergence of the integral of equation 



( ^.15 ) is automatic and the resulting field J-'(A) is everywhere smooth on space-time. 

The linear operator : X t-^ is defined on K[x] and we denote its range by T[x]. Also denote by F, r+ and 

r+[a:] the images under the operator !F of the spaces K, and respectively. 

It is clear that the space F[a;] consists of certain polynomials in the variable x°' with coefficients in the space F, so 
to understand the range of the operator T it is sufficient to identify the space F. 

To this end, we first introduce for any a and /3, solutions of the wave equation in spacetime the n — 1-form 
Lo{a,(3) = a{*dp) — P{*da), where * is the Hodge star operator on forms for the given Lorentzian metric. Since the 
wave equation for a field (p may be written d * {d(p) — 0, it is clear that the form u!{a, (3) is closed. Define r2(a, (3) to 
be the integral of the form w(a, /3), over a spacelike hypersurface, oriented towards the future, asymptotic to spacelike 
infinity, for given fields a and /3, which are required to be such that the integral converges and is independent of the 
choice of that hypersurface. Denote by W the space of all solutions of the scalar wave equation with initial data, 
on any spacelike hypersurface, asymptotic to spacelike infinity, in the Schwarz class (the initial data for a solution 
on a hypersurface is by definition the restriction of (p and ^dff) to that hypersurface). Denote by M[x\ the space of 
all polynomial solutions of the wave equation and by M'[x\ its dual space. Then for each 0(x) G W, we obtain an 
element ^(c/)) of the space M'[x\, defined by the formula fj,{(j)){f) = ri((/>, /), for each / e AI[x]. This gives a moment 
map /i : ^ M'[x], (f) i— > /i(0). Then we have the following result: 

Proposition V.4 F = Ker{fi). 

The proof of this result follows immediately from the Fourier inversion formula. 

Note that the information in the moment map /i is completely contained in the formal power series defined by 
the formula: p{4>){pa) = r2(e*^°^ , (p), where the exponential e'^"^ is understood as a formal power series in the null 
covector pa- The quantity p{4>){Pa) is then a formal power series whose coefficients are tracefree symmetric tensors, 
representing the various moments of the field (p. In terms of initial data, the quantity p{(j)) represents all moments 
of the data for the field (j). In this language the space F is the subspace of the space W consisting of all fields 0, for 
which p(0) = 0. 

Our final aim in this subsection is to determine the kernel of and to prove the 



Proposition V.5 For A e K[x] 



T{X){x) ^ J e''"'^''X{x,k)E = for all x (5.16) 
^ X{idk+tk) = 0. (5.17) 



Here the scalar t is an indeterminate. Also in writing equation ( 5.17 ) it is to be understood that in each term of the 
expression of the function X{x, k) as a polynomial in the variable a;, the expression is ordered by placing all the factors 

of the variable x to the left, before replacing the variable x by the operator idt + tk. 

Before we proceed to the proof of the proposition we want to clarify the structure of equation ( 5.17| ) with an 
example. In the case A £ -ftr2[2;], we have the expression X{x, k) = x°'x^ (3ab + x°'Pa + /3, for some symmetric tensor (3ab: 
vector pa and scalar /3, each of which depends only on the variable k. For this case equation (5.17) reads as follows: 



= {idl + tk''){ldl + tk'')pab + {idl + tk'')Pa + P- 



(5.18) 



Note that the commutator (id^ + tk°-) {idi. + tk^) — {id'^ + tk^){id'^ + tk"') vanishes identi cally , so there is no factor 
ordering problem. Expanding equation (5.18) in powers of the indeterminate i, equation (5.18) is equivalent to the 
following three equations: 
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= k'^k^Pab, 

= id^k' + k-dl)Pab - ik'^Pa. 
= dldlPab ~ idWa - p. 



(5.19) 
(5.20) 
(5.21) 



Note that equations (5.17) -(5.21) implicitly require that one extend the function X(x,k) off the null cone of momentum 
space before writing these equations since the formulation of the equations uses the full derivative operator dk- However 
it must be possible to rewrite the equations so that they are purely intrinsic to the null cone. In the case of equations 
(5.19) - (5.21), one can see easily that these equations are equivalent to the following three equations: 

k k /3ab: 

{L'^e + fcU'')/3,fc - Zp'^epab - ip'k,k^(3a, 








(5.22) 
(5.23) 



(5.24) 



+ tp^'kaP^Pb + {l/2)p''pa{lhl3'' + Pi) - [p^kafP. 

Here the operator L° is defined for any fixed vector p° by the relation: 

= 2pbfc[^9;!l. (5.25) 

It is clear from its definition that the operator L ° is i ntrinsi c to the null cone. Then equations ( 5.22| )-(5.24) hold 
for arbitrary vectors p°' . Alternatively equations (5.23) and (5.24) can be rewritten without using the vector p° as 
follows: 



{L, "A:^ + fc"L, '')l3ab - 3A:Vc6 - ik^k'^Pa, 



Pab 



ik(,La) "Pt + 2/3,rf + z/c(,/3rf) + {l/2)g,a{ikbP'' + Pb") - Kk^p. 

b] 



(5.26) 
(5.27) 



Here we have introduced the intrinsic operator L"'', given by the formula: L"'' = 2fc[°9j, , in terms of which one has 
the relation L°- = PbL'"^. Note that provided that equation (5.22) also holds, each of the equatio ns (5.26 ) a nd (5.27 ) 
amounts to just one scalar equation, since one may verify that the right hand side of equations (5.26) and (5.27) are 
proportional to the quanti ties k ^ and kck^, respectively. 

Proof of the proposition (5.7'/): We shall prove the statement for A € [a;], for every non- negative integer m 
by induction on the natural number m. First the required result holds for m = 0, since in this case the function 
X{x,k) — P(k) for some function P{k) G K. When equation ( 5.16| ) holds, the integral / e*'^"^ P{k)'B. gives the zero 
solution of the wave equation □$ = and it is well known in this case that this entails that the function P must 
vanish identically 

Next suppose the required result is true for a ll A G /^m[2;] for all m < s, for some positive integer s. We prove 
the result for m = s. So consider equation (5.16) with the function A(a;, A;) € Ks[x] now a polynomial in the variable 
X of degree not more than s. Then we can decompose the function A(a:, fc) as X{x,k) = a{x,k) + P{x,k), where 
a{x,k) e K<^'^[x] and P{x,k) e Ks-i[x]. 

Applying the wave operator to equation ( ^.16| ), we get the following equation: 



= J e'''-''\2ikad'' + n)\{x,k)E. 



(5.28) 



Since the function {2ikad^ + □)A(a;, k) E A's_i[a;], we get by the inductive hypothesis, the equation: 

= mkad^ + □) X{x, fc)]_(,a,+,,) . (5.29) 
Next take the partial derivative of equation ( ^.16 ) with respect to the variable a;. We get the following equation: 



= 



/ ikbX 



{ik" +d''){a + P)E. 



(5.30) 

Now we have a{x, k) = x'^ae{x, fc), where the function a'^(x, fc) = s^^d'^a(x, k) G Ars-i[x]. Then equation ( ^.30 ) may 
be rewritten, using an integration by parts as follows: 



ikhX 



{ik''x''ae+ik''P + d''X)E 



g«fc.x^(2zfc['^a;"lae + «fc''/3 + a''A)S H 



ix"'k^ae 5 



(5.31) 
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Now using equation ( 5.28| ), we have the following: 

Hence one has the following equation: 

J e''=''^'"ifc"aa5 = J e"'"''''' {2ikad''f3 + D\)E. 
This equation is used to replace the last integral of equation ( |5.3l| ). We then obtain 

= / e''=''^'"(-2fc['^a^lae + ik^P + 9" A) - |^(2ifc„c'"/3 + DA) S. 



(5.32) 



(5.33) 



(5.34) 



Now, by inspection, each term multiplying the quantity e*'^''^ of equation ( 5.34 ) is of degree at most s — 1 in the 
variable x. Therefore we may invoke the inductive hypothesis again to deduce 







-2/j["9^1ab + ik^fi + d^X ~ ^{2ikad''P + DA) 



2s 



2ik^''x'''^ab + ik^f] + d^X {2ikad''(3 + DA) 

2s 



Now equation ( 5.28| ), when written out gives the equation: 

[2^ka^'^(3 + ^Xl^^^,^^^,^ = [~2^ka^^al^^^,^^,,^ 
Substituting equation ( ^.36 ) into the last part of equation (5.35) gives the following equation: 



(5.35) 



(5.36) 



= 



2ik^''x''^ab + ik''l3 + d^X + ix'^ha'' 



(5.37) 



In the transition from the penultimate to the last line of equation (5.37), we have used the fact that the terms 
arising from the commutator of the operator of multiplication by fc" and the operator idk + tk exactly cancel 
the derivative term, the quantity i(9"A. That this is correct may be seen as follows. Consider the quantity 



[(fc--z9-)((xV)"/(fc))] 



{idk -\-tk) 



for n a non-negative integer, p a constant covector and for f £ K. We then 



have the following equation: 

[(r - ((a.V)'7(fc))]_(,a.+*.) - [(-V)"^V(fc) - *«P''(xV)"-V(fc)]_(,,,+,,) 

^ {{idi+tk')p,rk-f{k) - hP''(2;%)"^V(fc)]_(,a,+,fe) 

= r + tk')p,Y f{k) + ^np'^ + tfc^)pfa)""' f{k) - inp'^ [(^%)"" V(A;)] 

= fc"[(-'p^)"/(fc)].^(.a.+*.)- (5-38) 

Since any polynomial A € K[x] may be written as a finite linear combination of terms of the form {x^pi,)"f{k), we 
have the relation [(fc*" - id'')X]x^(^ig^+tk) = fc° \Mx>- >(id,^+ tk), for any A G K[x], as required. 

Finally we remove the factor ifc° from equation (5.37) giving the required result and the induction is complete. 



C. The general solution of □"'t/) = 

In this subsection we provide the general solution of the equation □"'■0 = in the space T[x]. As the construction 
will show this is equivalent to finding the general solution subject only to the condition that the zrm field D™"^^ lies 
in the space F (i.e. has zero moment map). We begin by proving that there is no loss of generality in restricting the 
domain of the Fourier transform operator ^ from the space K[x] to its subspace L[x]. More specifically, we have 
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Proposition V.6 For each X{x,k) €E K[x] there exists a ^{x,k) €! L[x] such that T{X) = ^(/i). 

Proof: For any A G K[x\, we have the decomposition, directly analogous to that of equation ( ^.12| ) above: 



A7- — 



00 

E 

s=0 



(-1)^x2*' r(i/ - 2r - s){v - 2r) 

22s+r r(i^-r + i)r(s + i) 



(5.39) 
(5.40) 



Here we have put x^ = Xqx" and v = (n — 2 + 2x'^9a)/2. From equation ( 5.40 ), by differentiation, it follows that 
the coefficients belong to the space L[x\. In particular it follows that the space K[x] is the sum of the space 
L[x\ with the module generated over the ring K[x] by the function x'^ . In view of this decomposition, to prove the 
required result it suffices to show that for any function a G /^m[2;] there exists a function /3 S Km+i[x\ such that 
JF(x^a) = ^(/3). Now we have the integration by parts identity, valid for any Ua G K[x]^ such that = 0: 







(5.41) 



Rephrasing equation ( 5.41 ) in terms of the operator we have J-{d'^ya + ix'^Da) = 0, valid for any i/a G K[x], such 
that yak"" — 0. In particular consider the case that ya = —ikQ^{koXa — taX^ki, + kaX^ti)a, where i° is a fixed unit 
vector and fco = k°'ta- Then it is clear that j/q G A'„i+i[x] and that ya satisfies the identity yak'^ — 0, so one has 
T{—d°'ya) = !F{ix'^ya)- But by contracting the vectors and iya we also have the relation ix"j/a = x'^a, which yields 
the relation !F{x'^a) = J^(/3), where f3 = —d'^ya- Since from its definition it is clear that the function (3 lies in the 
space iCm+i[a;], the proof is complete. 

In view of this result we henceforth assume without loss of generality that is defined on the space L[x]. 

The wave operator of spa ce-tim e, □ acts naturally on the space r[a:] and one has the relation, immediate from the 
definition of in equation (5.15|), valid for any A G L[x]: 



DT{X) = 2iT{D\). 



(5.42) 



We next wish to determine the kernel of the operator □ acting on the space r[a::]. By equation ( ^.16 ) above and using 
equation ( ^.42 ), we have □JF(A) = 0, for A G L[x] if and only if the function \{x, k) obeys the equation: 



= {DX){id + tk,k). 
Let the power series expansion of the function X Cz L[x] be given as follows: 



00 ^ 

X{x, k) = ^ ^ —x'^^x"'^ . . . x"'^ Xaia2...ar.{k)- 



(5.43) 



(5.44) 



In equation (5.44), each coefficient tensor Xaia2...ar belongs to the space K and is completely symmetric and tracefree. 
Writing out equation ( 5.43] ) in terms of this expansion gives the following system of equations, one for each positive 
integer q: 







00 ^ 



p=0 



(5.45) 



Note that there are no factor 



Here the quantity Ap.g is by definition Ap.g = i^d'^^ . . .d'^^k''^ . . .k''''Xai...a 61... 
ordering problems for the quantity Xp^q, since the tensor coef ficient s are all tracefree. Now expanding in powers of the 
indeterminate t we have the identity, derived from equation (5.44): 



X{id + tk,k) = J2 -Tl^pAk) 



p,q=0 



p\q\ 



(5.46) 



Comparing equations ( 5.46| ) and ( 5.45| ), we see that if we define //(x, k) = X{x, k) — X{id, k), then we have the relation: 

pL{id + tk,k). (5.47) 



17 



Note that /i G L[x], so from equation ( ^.47 ), one has ^ G ker(^). This gives the relation ^(A) = where 



X{id,k) e ioN- So we have shown that i f DT (X) = 0, then A = mod (ker(J^) + Lo[2;])- Conversely, if 



A e ker(jF) + io[a;], then it is clear from equation (5.42) that □JF(A) = 0. So we have proved the relation ker(njF) 
ker(jF) + ioN- Rephrasing we have proved the relation ker(n) = !F{Lo[x]). 
This generalizes immediately to our first main result 

Theorem V.7 For any integer m > and any space-time dimension n > 2, the kernel of the operator when 
acting on T[x] is given by the relation 

ker(n'"J^) = ker(J') + Lm-i[x]. (5.48) 

Equivalently this may he stated as: 

ker(n'") -^(L„_i[x]). (5.49) 



Proof: The required result has just been proved in the case m = 1, so henceforth assume to is a fixed integer 
greater than one. Suppose that $ G ker(n™) and $ = JF((/)), for some (j) G L[x]. Then D™^!*!) e ker(n), so 
we have = T[a), for some a G Lo[a;]. By proposition ( 5.73| ) the operator D is surjective as a linear map 



from the space Ls[x\ to the space Ls_i[x], for any positive integer s. It immediately follows that the operator 
is also surjective as a linear map from the space Ls[x\ to the space La-r[x\, for any positive integers r and s, 
with s not less than r. Therefore we may put a{x^k) = {2iD)™~^f3{x,k) for some /3 G Lm-i[x]. Then we have 
= = jc- ((2iD)™-i/3) = □™-i(J^(/3)). So we have $ - J^{(3) G ker(n'"~i) and the required result now 

follows immediately by induction. 

We have shown that the general solution of the equation □'"+^0 — in the space r[a;] is given by an integral 
formula 

<f>ix) = J e'"'^" I A(A:) + A,,...,, (fc)a:'^^ . . . x*^' | S (5.50) 

where the polynomial inside the integral satisfies the wave equation. The solutions of this type are automatically C°°. 
So with this formula we can not reach functions which are merely C*^ differentiable. However, this class of solutions 
is sufficient for our purposes. More general solutions can be obtained using functional analytic methods by starting 
with this integral formula on an appropriate function space and then taking limits. We will not pursue this here. 



D. Spinor momentum space 



We now specialize to the case of four dimensions and introduce spinor variables. A future pointing null momentum 
covector ka may be factorized as ka = k^kA', for fc^ a two component spinor, with complex conjugate spinor fc^/. 
More precisely, we have a surjective map from the momentum spin space to the future null cone of momentum space, 
which maps the spinor kA to the null vector kAkA' ■ The inverse image of the vector ka = kAkA' is the circle of spinors 
akA, (with a G C and |a| = 1) for ka non-zero, and is the zero spinor only, when ka = 0. We pull back our previous 
constructions along this surjection. The puUback of the ring is then the ring of functions f{kA, kA') which are 
everywhere smooth and vanish to all orders at the origin, decay faster than any power at infinity and which obey the 
differential equation {kAd^ — kA'd^ )f — 0. We shall need to multiply by components of the spinors kA and fcyi', so 
it is natural to enlarge the ring in the spinor case to the ring K which is by definition the ring of all functions 
f{kA,kA'), such that / has a decomposition as an infinite sum: / = X]r=-oo /f > with only a finite number of the 
functions fr non-zero and such that for each integer r we have: 

(i) fr is globally defined and smooth on the momentum spin space; 

(ii) limt^oo /r(e*"A;^, e*" kA') — 0, for all a G C, such that 3?(a) is non-zero; here the limit is taken with t real and 
the limit must be uniform on compact subsets of the momentum spin space; 

(iii) {kAd^ - kA'df)fr{kA, kA') = rfr{kA, kA')- 

For each integer s, denote by K'^ the subspace of K consisting of all / G if with fr vanishing for all r different from 
s. Then the puUback of the ring if^ is the ring K'^ and one has K^K'' included in Rp^'^, for all integers p and 
q. In particular K° is a subring of K and every space is a if^-module. Denote by K^x] the subspace of K[x] 
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consisting of all polynomials in x with coefficients in the space . Denote by L[x\ the subspace of k[x\ annihilated 
by the wave operator □ and by L^[x\ the intersection of the spaces i^-'[a:] and L[x\. Denote by if [tt] the space of all 
polynomials in the spinor variables it a and tta' with coefficients in the ring K. For p and q any non-negative integers 
and for j any integer, denote by Kp q[K\ and K-j^ ^[Tr] the spaces of all polynomials in the spinor variables tta a-nd t^a'^ 
homogeneous of degrees (p, g) in the pair (tt^jTt^'), with coefficients taken from the spaces K and , respectively. 
For every element f{x) of the space L[x\ and 5(71) of the space -?C[7r], we have unique expansions of the following form: 

/(^) = E/^i^-^^^i^i-^;^^'^'^^^'^^---^'^''^''-' (5.51) 

r=0 

00 

5W= E 5AiA....A,AiAi...A',^-^^7r^=...^^''7r^'^^^^...7r<. (5.52) 

p,q=0 



In equations (5.51) and (5.52), the coefficient spinors fAiA2...A^A{A'^...A'^ and gAiA2...A^A{A'^...A'^ lie in the space K 
and are completely symmetric in all indices. Denote by E'^ : L[x] — > KItt] the evaluation operator which substitutes 
the spinor tt^tt^ for x° in any element of L[x\. Then it is clear that the map is an isomorphism of the space 
L[x\, with range the subspace of KItt] consisting of all polynomials g(7r) S Klir], which obey the differential equation 
{nAd^-^A'd^')g{^)=Q. 

To proceed we need the spin or an alogues of our previous technical results. First the analogue of the surjectivity of 



the operator D of proposition (V.l 



Proposition V.8 The operators kAd'^, kA'd"" : L ^ L and kAd^, kA'd^ : K K are surjective. 

Proof: Using the isomorphism E^ it is easily seen that it is sufficent to prove surjectivity for the operators kAd^ 
and kA'd^ . Further by formal conjugation the proof of surjectivity for the operator kAd^ will yield a proof of 
surjectivity for the operator kA'd^ . So we just need to prove that the operator kAd^ is surjective, when acting on 



the space /^[Tr]. Using the expansion of equation (5.52), we reduce to proving that given a totally symmetric spinor 

gB...CB'...C' e K, there exists a totally symmetric spinor fAB...CB'...c' e K, such that k^ fAB...CB' ...c = 9b...cb' ...c ■ 
By taking components with a fixed primed spinor basis, we reduce further to the case that the spinor gB...CB' ...c has 
only unprimed indices. By contracting throughout with a spinor variable n^, we reduce to solving the differential 
equation kAd^f = g, given g e K[Tr], such that the solution / lies in the space K[tt] and both / and g are independent 
of the variable tta' ■ 

Let t° denote a fixed unit timelike vector and put ua = tak"^ . Note that UAk"^ = tak^k^ is always a positive real 
number unless kA — Q- Then one has the following decomposition of the function g: 



^-^{kATT^nnBTT''ygp.,q. (5.53) 



^-^ plql 

p.,q=0 ' ^ 

This decomposition follows from the expression of the spinor t:a in terms of the spinor basis ua and fc^: t:a — 
{f^kckc')^^ {—ksT^^ riA + nsTr^/c^). Then by the binomial theorem, we have the following explicit formula for the 
quantities gp,q: 

^ [t^kcLy^^ [^^Ad^rikBd^yaU, ■ (5-54) 

It is clear from equation ( ^.54| ) that each coefficient gp^q lies in the ring K, so by linearity it suffices to prove the 
required result for the case g — (fc^7r"^)^'(nB7r^)'?, with p and q non- negative integers. But then we have the following 
relation: 



kAd^ 



{kA^^Y{nB^''Y+^ 
{q + Ijf'kckc' 



= ikATT^rinBTT^'y. (5.55) 



So / = {{q + Vjf^kckci) [kAi^'^Y {nBi^^Y^^ provides a solution in this case. Since it is clear that this function 
/ belongs to the space K[k\ and is independent of the variable tta', the proof is complete. Note that by tracking 
homogeneities through the proof we find that if g belongs to the space if^ then we may take the solution / to lie 

in the space K^'^l^^. 
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Second we need to analyze the kernel of the pullback of the Fourier transform operator T . This Fourier transform, 
still called JT, is defined now as follows, when acting on any (l> ^ K\x\: 

T{(j))= [ e'''-''°(j){x'',kA,kA')n. (5.56) 



Here one has f2 = g-^^g^ ^ dkAdkBdkA'dks' and the integral is carried out over all of spin space. It is easily shown 
that the operator !F maps K'^ isomorphically onto r+ (the range of T of section 5.2 acting on the space K~^) and 
annihilates all the spaces , for j non-zero. Furthermore, acting on the space i<r°[a;] the operator !F agrees with 
the pullback of our original Fourier transform operator (restricted to the domain i4r"'"[a:;]), up to a fixed non-zero 
multiplicative constant. 

Proposition V.9 ; For each cj) e L°[a;]; 

(PekciT ^ (P{x'',kA,kA') = {d^ + ix''kA')<l>A + {df +tx''kA)<l>A', (5.57) 

with (f>A € L^[x\ and (pA' G obeying the spinor zrm-field equations: d°'4>A — and d'^cfiA' =0. If (j> has degree 

at most m in x, then (/ja o,nd (f>A' rnay be taken to have degree at most m — I in x. 

Proof: The "if" -part of t his re sult is a trivial integration by parts, so we assume that J-^cj)) vanishes and we establish 
the formula of equation ( ^.57 ) for the function (p. First if ip is independent of the variable x, then J^((/)) = entails 



that (f> = 0, so the result holds if we take (pA — (pA' — 0. So now we assume that the required result is true for (f> any 
polynomial of degree at most to — 1 and take to have degree at most to. Then applying the wave operator to the 
equation J-{(p) = 0, we obtain the equation T{kad°'(p) = 0, so by the inductive assumption, we have the relation: 

= {d^ + ix''kA')tpA + {df + ix''kA)i'A' ■ (5.58) 

Here the quantities ipA G ^^[x] and %pA' G -^^^H are polynomials in x of degree at most m — 2, belong to the spaces 
L^[x] and and obey the field equations d"''4}A — and d""ipA' — 0, respectively. Write (p = a + (3, where a is 

homogeneous of degree exactly m and [3 is of degree at most m — 1 in the variable x. Similarly decompose the fields 
ipA and ipA' as ipA = —ipA + cta and -ipA' = —ipA' + ca' , where pA and pA' are homogeneous of degree m — 2 , whereas 
a A and aA' have degree at most to — 3. Then the terms of highest degree in the variable x of equation ( 5.58| ) give the 
following equation: 

kad^a = x^kA'PA + x'^kAPA' ■ (5.59) 

Note that by the inductive hypothesis, the functions pA €E L^[x] and pA' G "'^M obey the zrm field equations: 
d'^PA — and d°'pA' — 0. Note that the quantities a, p and p' are respectively of homogeneity (to, to), (to — 1, to — 2) 
and (to — 2, TO — 1) in the variables tta and tta', respectively. Now put a;° = tt^-k^ in equation ( p. 59 ). We obtain the 
equation: 

kad^d^'a = m{kA'Ti^'p + kA^^p'). (5.60) 

Here we have put a = E^{a) G if^,„[7r], p = E^in^pA) G i^™_i,™_2 W and p' = E^{n^'pA') G i^^La.m-iW- Next 
write p = m~^kAd^T and p' = m~^kA'd^ r', for some r e K^n.m-2[''^] ^^'^ G -f''m-2,mM- This we can do by the 
surjectivity of the operator kAd^ proved above. Then equation ( 5.60| ) may be rewritten as follows: 



= k,[d^d^'a - TT^'d^T - TT^'dfr'). (5.61) 

Now suppose that the quantity v°- e K[tt] obeys the equation kaV°- — 0. We may expand the vector in terms of 
the spinors k^, and their conjugates fc"^ and n"^ as follows: 

v'' ^ k^k^'U + k^n^'V + n^k^'w + n^n^' X, (5.62) 

Ov^UAUA^ {v'^nAkA') , . 

{t.kCkC'Y' ^- (tjz^kC'Y' ^^-^^^ 

Jv^-kAUA.) {v^kAkA± , . 

{t.kCkC'Y' - {t^kCk^'Y ^ ' 
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It is clear that each of the quantities U, V, W and X lies in the space K[Tr]. When we have the relation v°'ka = 0, this 
implies that the quantity X vanishes. This in turn entails that the quantity w° may be expressed as = k^v^ v^, 
for some G K[tt] and G K[n]: indeed one may take = Uk^ + Wn^ and = Vn^ . Note that if 
v"^ £ K-^ ^[k], then by equation ( 5.64 ) the quantities and may be taken to lie in the spaces -/^^^^[Tr] and 
if^^^[7r], respectively. 



Applying this result to equation (5.61), we obtain 

d^dt'cj = TT^'d^T + ir^dfr' - m'k^v^' - m^k^'v^, 



(5.65) 



for some € „_]^[7r] and G i^mi'Li „i_i[7r]. Contracting equation ( 5.65 ) through with the spinors tta and 

TTA' gives the following equation: 



(7 — k^TTAv' + k^' TlA'V- 



(5.66) 



Here we have put v = ttav"^ G m-iH ^'^^ v' = tta'V"^ G K^^^_^ mM- Rewriting equation (5.66) in terms of the 



variable x, we find: 



a = x°-{kAaA' + kA'aA)- 



(5.67) 



Here the fields aA{x) and aA'{x) are determined by the formulas E^{a ) — 'm~'^d^v and E'^{a ) = to~ 9^ v' . Also 
we have aA^x) G L^[x] and aA'{x) G L~^[a;] and both the spinor fields a^4ja;] and Q!y!i'[a:] obey the zrm field equations 
and are homogeneous of degree m — 1 in the variable x. By equation ( ^.67 ), we have the following relation, using an 
integration by parts: 



= Tic/)) = Tia + I3)= Tix^ikAaA' + kA'UA) + P) 
= J'iidfaA' + id^UA + 13) 
By the inductive hypothesis, we obtain from equation ( 5.69| ) the relation: 

idfuA' + id^aA + (3 = {d^ + ix^kA'pA + {df + ix°-kA)uJA' 



(5.68) 
(5.69) 

(5.70) 



Here the fields oja G L^[x] and lua' G i ^[x] obey t he sp inor z rm fi eld equations and are polynomials of degree at 
most m — 2 in the variable x. Combining equations (5.67) and (5.70), we get: 



= a + /3 = x^'ikAaA' + kA'aA) + [d^ + ix''kA')oJA + (9^' + ix''kA)LOA' - i{dk a.A' + d^aA) 



{d^ + ix''kA')(t)A + {df + ix''kA)c^ 



(5.71) 
(5.72) 



In equation ( 5.72 ) we have put 4>a = — iotA and (pA' = uja' — ictA- Since it is clear that the fields ipA and (/)a' 
have all the requisite properties, we have proved the validity of equation (5.57) for any field (t>{x) G L^[x] of degree at 
most TO in the variable x. Therefore by induction we have the validity of equation (5.57) in general and the proof is 
complete. 



E. The general solution of M'$ — 

If we wish to construct a space-time field from elements of with j non-zero, we first need to multiply by spinors 
kA or kA' as appropriate to map the element to an (indexed) element of K^, before applying the Fourier transform 
operator J-'. The result is a spinor indexed field on space-time. For example consider the standard zrm equation 
qAA ^j^g (jjji^x) — for a totally symmetric spinor field ^ab...cd{x) of r indices. Taking another derivative and 
contracting, we imm ediately find that the field ^ab...cd{x) obeys the wave equation D^ab...cd{x) = 0. Therefore 



by the theorem ( V.7 ), we may write its general solution, (after pulling back to the momentum spin space) in the space 
r+ , as follows: 

'i>AB...CD{x) ^ j e'''-'''^aAB...CDikE,kE')^. (5.73) 

Here the Fourier coeflicients aAB...CD lie in the space K*^. Applying the field equation we get the equation 
k^aAB...CD = 0, whenc e it fo llows that aAB...CD{kE, ks') = kAkB ■ ■ ■ kckD4>{kE, ks'), for some function (j){kE, kE') G 



K ^ . So now equation (5.73) reads as follows 
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^AB. 



CD 



ix) 



kAks ■ ■ ■ kckD4>ikE,kE') ^■ 



(5.74) 



Next we shall derive explicitly the solution by Fourier transform of the equation d^'^^ ^ab '^c = in the special case 
of one primed index (m = 1) and then later ge neral ize to arbitrary p ositi ve m. We know that the field ^abb' lies 
in the kernel of the operator by proposition (H.l), so by theorem ( |V.7| ) it admits a Fourier representation of the 
following form: 



$(X) 



(abX^ + a) 



(5.75) 



Here the variable X is an abbreviation: X = (x", tt^, tt^') and we have put $(X) = ^ABB'{x)TT^Tr^TT^ . The Fourier 
coefficients and a depend on the spinors fc^ and tta and their conjugates, but not on the variable x"'. Defining 
the operator M' = tt^ d^da the field equation may be written as M'$ = 0. This operator agrees with the operator 
M' defined in section 2 in its action on the spinor indexed coefficients of $. Applying the field equation, we get the 
following equation: 



= T (^TT^'af (ifcb + db){acx'' + a)) =T {(t:^' df){iki,acx'' + ik^a + ai,] 



Using equation (5.5?) above, we deduce the following equation from equation ( 5.76| ): 

{n^' d^){ikhacx'' + ikta + at) = {d^ + ix°-kA')aA + {df + ix°-kA)aA' 



(5.76) 



(5.77) 



Here the quantities a a and ua' are independent of the variable x. Equating the coefficients of x in equation ( ^.77 ) 
we get 



i-K^ ksd^ac = kcac + kcac- 



(5.78) 



Equation ( 5.78 ) gives immediately the equation ksdil^a'^kc ~ 0, which is solved by a'^kc — {tt ks) oe, for some a, 



independent of the variable tta- This gives the relation: Oc = fc^Tr ttc/^ c + kgic + kc'^c, for some f3c', Ic and 



7c, with a = k^ (3a'- After an integration by parts applied to equation (5.75), the terms involving ^c' and 7c may 



be eliminate d, so one may take just ac = ksTT^TTcPc' , without loss of generality. 
Equation (|j^ now becomes 



{n^' d^){ikbkDTT^ncPc'x'' + iha + koTT^TTEPB') 
= iks'T^^'koT^^kcPcx" + ikB'TT^'kgd^a + Skon^TT^' f3B' 
= {d^ + ix''kA')aA + {df + ix''kA)aA'- 



(5.79) 



This gives the relation fc^/a^+fcyiaA' — ks'TT^ kuTT^kAPA'- So we have a^' = ks'TT^ knir^ ^A'+^kA' anda^ — -_kA5 
for some 5. But then the contribution of the terms involving the quantity S to the right hand side of equation ( 5.79| ) 



is just {d^kA — d^'kA')5 = [kAd^ — kA'd^ )S, which vanishes, since by tracking homogeneities we find tha t S li es in 
the space if". Therefore without loss of generality, we may take 5 = 0. Then if we put 
obtain the relation: 



^A' 



in equation (5.79), we 



ikE'Tr"" kBd^a + 2kDTT^TT^' Pb' = fes'^r^ fcuTr^Sf /3a 



B'l 



^DciA' , 



(5.80) 



Next write (3a' = Pa'B'T^^ + t^A'P, where Pa'B' is symmetric and both (3a'B' and /3 are independent of the variables 
TTA and TTA'- Then putting tta' = ^a' in equation (5.80) gives the relation: ^A'B'k^ k^ = 0, which entails that 
Pa'B' — kA'^B', for some spinor Sb'- Then /3a' = kA'^B'T^^ + tta'C, for some scala r e. The term in Pa' proportional 
to kA' may be eliminated by an integration by parts, applied to equation (5.75), so we may take without loss of 
generality: (3a' = tt^'C. Then equation (5.8C) reduces to the equation: 



ikBd^a = koTT^TTA'd^ e 



(5.81) 



Next we put a = OAsTr'^Tr^, for some symmetric spinor uab, independent of the spinor t:a- Putting tta = kA in 
equation (5.81) gives the relation: aAB^^k^ = 0, so a^s = kA^B: for some spinor e^. Also put eb = eBB'^^^ , where 



eb is independent of the spinors tta and tta'- Then equation (5.81) reduces to the following equation 
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(5.82) 



Then we have the relation a^x^ + a = k^Tr^ {x^ttbttb'£ + ^bb'T^^t^^ )■ Summarizing we have found that the general 
solution by Fourier transform of the equation d^^^ ^ab = is given by the formula: 



$(X)= / e**=°^"fcB7r^(x%c7rc'e + ecc'7r^7r'=")J7. 



(5.83) 



Here the quantities e and ^bb' depend only on the momentum spinors kA and kA' and are subject to equation ( ^.82 ) 



Next equation ( 5.82] ) may be solved by writing e = — z0^fc"^, for some (pA- Then one has e^"^ = i9j 



A' aA _ i.AaA^ 



for 

some (j)"^ . The quantities (pB and (j)^ are freely specifiable. Writing out equation (5.83) in terms of the spinors (pA 
and (jjA' J we get the following formula: 



$(X) = J e''''"^°fcs7r^(-ix%c7rc'0B^^ + (^f^^ - k^ (f)^' )ttct^c')^- 



(5.84) 



Finally we may integrate by parts in equation (5.84) to eliminate the derivative term. This gives the equation: 
$(X) = j e'-^'^'^'' kBTT^ {-ix'-^'^' TTcT^cxpBk^ + ix'-^''^ kc4>B'^^ t^c - k^ nccj)'^" ttc) ^ 



{kBn^)^{ix^'^(f>cnc' + TTc) ^ 



c' 



(5.85) 



Note that there is gauge freedom in the pair of Fourier coefficients (pa 
is unchanged under the transformation < 

0)A' 



^ ): the quantity e of equation ( 5.82| ) 



-> (pA + fcyiT- Then the quantity e is also unchanged provided we make 
the transformation (p"^ i— > cp"^ + 7. So the complete gauge transformation is (pa ^ (pa + Kaj, where Ka is the 



operator pair: Ka = ( kA, ). 



Note that equation (5.85) may be rewritten in the following compact form: 

^{X)= j e'^'^''\kBTT^fZ"(Pan. 

Here Z" is the twistor Z" = {ix""iT a' t A') ■ This result may be generalized immediately: 
Theorem V.IO The integral 

$(X)= /e*'=»^>(Z",/cB7r^)0, 



(5.86) 



(5.87) 



is a solution of the equation A/'$ — for every (p{Z°' , kBT^^) G K[Z, kATT^]- Conversely, let $ G r+[x,7r,7r] be a 
polynomial homogeneous of degree {p,q) in the spinors (7r,7f) with p > q. If ^ satisfies the equation M'$ = then 
it has the representation (5.81) for some (p{Z, kAT^^) G K[Z, kAT^"^] homogeneous of degree {q,p) in the variables 

{Z,kATT^). 

The restriction p > q is necessary, because the theorem is false when p = q, ot when p < q, because in each of these 
cases the equation M'<I> = possesses a gauge freedom: if $ = {■KA'''^Ad°'YP: where p is a polynomial homogeneous 
of degree {Q,q) in the variables [ttatT^A')! then the equation ir^ d^da^ = is automatically obeyed, for arbitrary 
such functions p, as is checked easily , since the operators Tr^'Tr^c?" and M' commute and since p is annihilated by 
the operator M' . Because of this gauge freedom, no Fourier transform formula based on the null cone of momentum 
space is possible, unless one first fixes the gauge freedom in some way. 

Proof: By straightforward differentiation we see immediately that the function $ obeys the equation tt^ d^da^ = 0. 

Conversely we prove next that the general solution of the equation Af'$ = may be put in the form of equation 
(5.87). The proof is by induction on the integer q. First consider the case g = 0. Then the field ^{X) is independent 
of the variable tta' i so the field equation M'$ = is equivalent to the equation d^da^ = 0, which is just the standard 
zrm field equation, for a totally symmetric spinor field with p > Q indices. By equation (5.74) above the solution may 
be written as follows: 



^X)^T((kB^' 



(5.88) 
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Here the function (f) is independent of the variables {x°' , it a, t^A')- Therefore the required result holds in this case, with 
the function /{Z") independent of the variable Z". 

Next consider the case q > 0. Let ^{X) of homogeneity {p,q) satisfy the field equation — and put 

^'(X) = d:^ d^da^- Then ^{X) is of homogeneity (p — 1, g — 1) so, since q — 1 is non-negative and since p— 1 > g — 1, 
we may use the inductive hypothesis and write the field '^{X) as follows: ^' — J-{ikaX"'{kBTT^)P~^'ip{Z°')). Define a 
field F{X) by the following formula: 

(5.89) 

We wish to choose the function /(Z"), such that the field F{X) is of homogeneity {p,q) and obeys the equation: 
"if — d^daF. Applying the differential operator d^da to equation ( 5.89| ), we get 



d^'d^daF{X) ^ Hidfd^ida + lkMkB^''Yf{Z'^)) 
= Tid^'d^daikB-^^YfiZ")) = -pT ([kB-K^f-'d^'k^dafiZ")] 
= lJ'{{kBTT^r-'dfnA'kA{dz)^f{Z'')) = ip{q + l)J^{{kBTT^r-^kA{dz)^f{Z^)). (5.90) 

Therefore we have a solution to the equation {d^ d^da)F ~ 5* provided that the function f{Z°') obeys the equation 

kAidz^fiZ'-) = iip{q + l))-\g(^"). 

Having found the function F{X), we write ^{X) — F{X) + g{X), for some function g{X). Then the function 

g{X) must obey both the equations {d^ d^da)g = and tt^ d^dag — 0. Combining these equations, we get the 
single equation: d^dag = 0. By contracting this equation on the left with the operator ttbO^^ , we see that the 
function g{X) lies in the kernel of the operator □ whence it admits a representation analogous to that of equation 



(5.74): g{X) = ^^(7), for some function 7(fcA, fc^/, tts, tt^/). Then the field equation d^dag — gives the equation 
kA^^l = 0, so 7 = (fc^7r'^)P?7, for some function 'q{kA, fc^sTr^')- Putting (t){Z°') = f{Z") + 1] now gives the desired 
representation of the function $ and the complete proof follows by induction. 

If we apply this theorem to our special case we obtain the explicit representation given in the following 

Corollary V.ll For any non-negative integer m the integral 

i^B^.C' {x) = j e''="^"/c-^/c^ ...k^ {0B'...C' + <I^E(,B'..x''c') ' ' ' + ^I^E ...Fx"" B' ' ' ' c ] ^, (5.91) 

is a positive frequency solution of the equation dA(A'4'B^"ci) ~ ^' Conversely, every such solution is represented in 
the above form. 



VI. TWISTOR SOLUTION OF THE FIELD EQUATIONS 

Let us first introduce the structures which are relevant for the purposes of this work. For more details see |l^] 
and references therein. Twistor space T is by definition a four dimensional complex affine space. Denote by V the 
underlying vector space of T. At any z S T, denote by 9{z) the natural isomorphism of the tangent space of T at z 
with the vector space V and denote by 9 the V- valued one form on T whose value at any z S T is 9{z). It is clear 
that the one form is exact: 9 — d(. Here the quantity C is a V valued function globally defined on T. The function 
C is unique up to the transformation: C i— > + a, with a constant. The function (. serves as a vector valued global 
coordinate for the space T. 

Naturally associated to the affine space T is the space S{T) which is the space of all two dimensional affine subspaces 
of the space T. Naturally associated to the vector space V is the space M(T) which is the space of all two dimensional 
subspaces of the space V. There is a natural surjection, p : S{T) M(T), which takes each element of S{T) to its 
tangent space. The map fi renders S{T) a two dimensional fibre bundle over M{T) with fibre a two dimensional affine 
space. The space M(T) is provided with a natural holomorphic conformal structure, which is such that x,y G Af(T) 
are null related if and only if the intersection x H y is non-trivial. It is isomorphic to (compactified, complexified) 
Minkowski space and the space S (T) can be regarded in a natural way as the affine (unprimed) spin bundle over the 
space M(T). The primed cospin bundle S"(T) is by definition the space of all pairs (X, z) G S{T) x V with z tangent 
to X. It is a two dimensional vector bundle over the space S{T). Denote by L'{T) the line bundle 17^(5" (T)) over 
S{T). Note that the restriction of the form 9^ to any X e S{T) naturally takes values in the line bundle L'(T) at X, 
pulled back to the space X. 
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Let / denote a holomorphic function defined on some domain U in T. Then the form fO"^ is a holomorphic two 
form on T with values in f2^(V). For suitable X G S'(T), consider the following contour integral: 



s{f){x) = / fe\ (6.1) 

Here ^{X) is a closed oriented contour of two real dimensions, which is required to lie in the intersection of the 
space X with the domain of the function / and to vary smoothly with X. It is clear that the quantity S{f) represents 
a holomorphic section of the line bundle _L'(T) over its domain of definition, M{U) (this domain is an open subset of 
the space of subspaces X, for which the intersection with the open subset U has non-trivial second homology). By 
definition, if the integration contours are regarded as given, the section S{f) is the (unprimed) spinor field associated 
to the twistor function /. 

More generally let F{z, C) denote a holomorphic function on the tangent bundle of the domain U with (z, () E U xY. 
Then consider the following contour integral for {X, z) E S"(T), such that X e M{U): 

S{F){X,z)= [ F{zX)e\ (6.2) 
This defines a function S{F) on S"(T) taking values in the line bundle i'(T) (pulled back to the space S"(T)). 



Lemma VI. 1 The contour integrals of equation (6.1) and {6.i) give coordinate independent formulations of solutions 



of the zrm equation and the equation Af'$(a;'^,7r , tt^i/) = tt 9^9a$(a;'',7r , tt^/) — 0, respectively. 

Proof: We use lower case greek indices for tensors based on the vector space V and introduce the standard represen- 
tation = (C"^, Qai) of a twistor C, in terms of an unprimed spinor C,^ and a primed cospinor Q^i . A point X, not at 
infinity, of the space SiT) is labelled by the pair [x"" , tt^). The two dimensional afhne subspace of T corresponding to 
X is the set of all twistors C" {pA' ) of the form : 

C{PA') = {ix^^'pA' + TT^, PA') - x"^'pA' + n", (6.3) 

for arbitrary cospinors pA'. Defining X°'^' = {ix"^^' , 5%), = (tt-^, 0) and X"'^ = X"'^'x'^^'ec'D' , we note that 
the restriction of the one form 0" to the space X is given by the formula: X*{9°') = X"'^ dpB' and therefore the 
restriction of the form 6'"6"3 to the space X is just X*{e"9^) = X°'l^d'^p, with d^p = {l/2)e^' ^' dp a' dp w ■ Therefore, 
the field S{f){X) of equation (|I|) factorizes as 5(/)(a;'=, tt^')"^ = X"^0(/)(a;'=, tt^ ), where we have the following 
explicit formulas for the function (j){f){x°' ,t:^): 

0(/)(x^ TT^) = / /(x"^ + n")dV 

Jp{X) 

f{tx^'''pB'+7r^,PA')d^p. (6.4) 

P{X) 

In equation (p^) the two dimensional contour p{X) lies in the primed spin space of the variable pA' and varies 



smoothly with X, avoiding the singularities of the integrand. Differentiating equation (6.4), we immediately obtain 
the zrm field equations in the form d^da(j){f) — 0. 



Similarly, the integral of equation ( |6.2| ) gives rise to a field '^{F){x°- ^-ka') given by the following formula: 

^{F){x\n^,^A')= I ^^(x"^'^B';^"^PB'+n")dV 

Jp{X) 

F{ix^'''TrB',nA'-:ix^'''pB' +n^,PA')d^P- (6.5) 

Differentiation of equation ( |6.5|) immediately gives the field equation AI'^{F) — 0, as required. 

Note that, depending on the properties of the twistor functions / or F, the fields 4>{f) and ^(F) may contain many 
different helicities or irreducible spinor parts. 

Denote by 0{p,q) the sheaf of germs of holomorphic sections of rank p totally symmetric covariant tensors on 
projective twistor space -P(V), taking values in the sheaf 0{q) (the sheaf of germs of holomorphic functions h{() 
homogeneous of degree q in the variable C). Such a section is described non-projectively by a tensor with p indices: 
faia2...ap{(")d"^ ®0°''^ (g) . . . (g) 6'"" , such that faia2...ap is Completely symmetric, holomorphic, homogeneous of integral 
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degree q — p in the variable C" and such that = C"VQia2 - ap(C")- Our main result is that the sheaf cohomology 
group H^(U, 0(j), 2p— 1)), for p a positive integer describes the general analytic solution to our higher spin equations 
(the spin is p+1/2), for suitable domains U in twistor space. This restriction to analytic solutions is not mandatory: 
by replacing ordinary cohomology by C.R. cohomology one can obtain non-analytic solutions from the twistor theory. 
We do not discuss this further here. 



A. Twistor description for the spin (3/2) case 

As before we begin with the spin (3/2) case and treat the general case later. In this case the object of study is the 
group H^{U,0{1, 1)). We shall use the contour integral description to get at the results. Each calculation that we 
do then corresponds, according to well established procedures, to an appropriate calculation using sheaf cohomology 
as described for example in the books of Penrose and Rindler. For H^{U, 0(1, 1)) we use functions /q(C), which are 
homogeneous of degree zero in the variable C. For a function g(C), homogeneous of degree zero, the corresponding 
spacetimc field gsc is spin one and is given by a contour integral according to the standard formula of Hughston: 



9Bc{x'')= dedcgiix''" pB-,PA')p'- dpc. (6.6) 

Here the operator Ob denotes the partial derivative with respect to the unprimed spinor part C,^ of the twistor variable 
C^. Also the one dimensional contour 7(2;) is closed, avoids the singularities of the integrand and varies smoothly 
with the point x. 

Next we need the explicit action of the twistor operator on the field gBC- Multiplication of 5 by C gives a 
function homogeneous of degree one for which the corresponding field C^g is spin (3/2) and is given by the following 
formulas: 

{C9)bCd{^'')^ I dBdcdDiC9)ilx''^'pF',PE')p'''dpG' 

= f [U^dcdDgilx"^' PF',PE')+X^''' pB'dBdcdDg{lx''^' PF',PE')] p'^'dpG' 
J-yix) ^ ^ 

= 3S%gcD - iX'^^'dB'BgcD- (6.7) 

In equation (|6.7| ), the twistor 6^ = ((5b"^,0). Also we have used the fact that inside the twistor integral the operator 
dh is represented by t he ope rato r ipb'Ob- Applying these results to the indexed function ffsiC), allowing for the extra 



index, equations (6.6) and (6.7) become 



<t>pcD{x'')= I dcdDfp{tx^'''pB',PA')pc'dp^'. (6.8) 

{C(t>)%CDE = i^i^liDE " iX'-''' dc'c4>pDE. (6.9) 



If we now impose the condition C"/q — 0, we see that the trace over the twistor indices of equation (6.9) must vanish 



From the right hand side of equation (6.9) this gives the following condition 



= 30CDB - iX°'^ dcC"l)aDE = 3(f>CDE + X^^ dcC'<l>ADE ~ idA'd^DE 

= S4'CDE ~ '2.4)CDE + dcC'{x^^ (pADE ~ T-<i>iiE) = "^CDE ~ ^ecioi^E) ~ idA'C'4'DE^ (6.10) 

with (paBC = {4>ABC, <t>DE)^ 4>ABC = if^ABC + f^AiBi^c) and tp^^ = (/)^^ + ix^ ^(pcDE, where the spinors (ji^^^, iI^abc 
and iI'de completely symmetric. 



Equation ( 6.1C ) gives in particular the equations ipcDE = '^'^A'(c4'de) ^^'^ '^A = (— i/3)5 ijjABB'- Hence the field 
4>ABC is completely determined given the field V'ab- Once the field (j)ABC is known, the field (f>^B '^an be recovered 
from the formula: 4>ABi^'^) = ~ '''^^ ^ 4'CAb{x'^)- So knowledge of the single field V'ab (and its derivatives) is 

completely equivalent to knowledge of the original field 4>aBC- Finally the field equation for the field (f^aBC is just the 



standard zrm field equation d (j^aBC — 0. This equation immediately implies (by straightforward differentiation) 
the equation d^'^^ V'bc ~ ^ ^^'^ conversely it is seen easily that the equation d^''^ i^BC ^ iniplies the field equation 



^<PaBC = 0. So we have established that the cohomology group H^{U,0{1, 1)), for appropriate domains U in 
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projective twistor space is isomorphic to the space of solutions of the equation d^^^ i-'^c ~ "-"^ corresponding 



domain in space-time, with the field V'sc being totally symmetric. Finally from the definition of the field '4'ab ^'^^ 



A' 



from equation (6.8), we have the following contour integral expression for the field ^^g'- 



^ab{^'') 



dAdBf" (ix^'^'pB' ,PA')+ ix'' ''dAdBfc{ix''" pb' , PA')pc'dp 



AB' , 



(6.11) 



7(x) 



Contracting equation ( |6.11 ) through with tt^tt^tta', and using Cauc hy's integral formula to reduce the integral to a 
one dimensional integral, we find complete agreement with equation (6.5), where the function F(z, C) = z^fa and the 
field ^{F){x"' ,TT^ ,TTA') is then a fixed constant multiple of the field t(j^g{x'')'K"^TT^ tta' ■ So we have the shown the 



Proposition VI. 2 There exists an isomorphism between the sheaf cohomology group H^(U,0{1, 1)) and the space of 
holomorphic solutions of the equation dAA'4'B^ = 0- 



B. The general spin case 

Consider the interpretation of the twistor cohomology group H^{U,0{p,q)). We assume that the integers p and 
s = q — p-\-2 are positive (later we shall restrict further by requiring that s > p). As discussed above, a representative 
element is an indexed function faia2...a {C)^ which is completely symmetric, holomorphic, and homogeneous of degree 
q — p \n the twistor variable C and such that = C"^/QiQ2- -a (0- The unprimed spinor field corresponding to the 
fimction faia2---a may be given as follows: 



Z",a;",7r^) EE / i^(Z", + n")^^^. 

JliX) 



3.12) 



Here we have put F(Z,C) = Z"^ Z""-" . . . Z"" f„^„^,,,a^{C)- The field (/-(Z", a;°, tt^) is a homogenous po lynom ial of 
degree (p, s) in the pair of variables (Z", tt'^), with coefficients obeying the zrm field equation by lemma ( VI.l). Our 
first objective is to obtain a formula for the action of the twistor variable (^^ on such a field. Denote the result of this 
action by {C,(j))^ ■ Then for this field we have the expression: [C,cf))^{Z, x, n) = X^'^ (pc' {Z, x, tt) + Il^(f>, where the field 
(j)B' is given by the following formula: 



cj)c'{Z,x,7r)^ [ pc'FiZ^.X'^''' pB'+T\^)dPp. 



3.13) 



Multiplying both sides of equation (3.12) by s + 1, we manipulate equation (6.13) as follows: 



{s + l)(Pc'{Z,x,^)=^Ad^ f pc'FiZ^.X^''' pB' +n^)(fp 

= 7r^ / pc"{{dc)AF){Z",X"''' pB' + U")d^p = -t7r^dAC'^{Z,x,7r). 

Jl(X) 



Summarizing we have established the formula 

{Ccj^f^^^ (-^[s + l)-^X^^'^a + 5^A) ^■ 



(6.14) 



(6.15) 



Next, we consider t he co ndition = C"Vcna2 .-ctp (C): which is equivalent to the condition — C°'{dz)aF. When 
applied to equation (3.K), with the field (j) replaced by the field {dz)a4' this implies 







(iX'^^'da -{s + 1)5^) {dz)p<l> = -^^x^^'9„(9z)b0 -{s + 1)^^(9z)b0 + iii^dad^' cj,. (6.16) 



Next we perform a change of variables and write the twistor ~ {Z , Za') as follows: — X" zb' + z ^bj so 
we have Z^ — ix°-ZA' + z^ and Za' — za'- Then the field (f> becomes a function of the variables {x°',tt^, z"), where 
z" = {z^, ZA')- Under this change of variables, we make the derivative replacements: da da~ izA'{dz)A, d^, 
{dz)A ^ {dz)A and — ix°-{dz)A- Equation (3.16) then becomes 
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ABA' 



{da - iza'Oa) dB4> - (s + l)7r^9B0 + iT^^ {da ~ izA'dA) (d^ 



„BA' ', 



= TT^'dAi-s + 1 + 7')0 + iT^'^dad^' (j}. 



(6.17) 



Here we have put Oa = = {dz)"^ and 7' = za'O"^ . Also define 7 = z'^Qyi. Note that the field obeys the 

equation (7 + 7^ — = 0. 

Equation ( |6.17| ) may be regarded as giving a partial propagation of the field (p in the directions. But we also 
have the field equation obeyed by 0, which in terms of the original variables is the zrm equation d^da(j) — 0. In terms 
of the new variables this equation becomes the equation: 



0^d^da<t>-tZA'd^dA(l>. 



Removing the factor ir^ from equation (6.17), we get the following equation, valid for some field 

dA{-S + 1 + l')P + ldad^'<i> ^ ^AX- 

Applying the operator to ( 6.19| ), we get 



{s + l)x = d^dA{-s + l + l') 
Then, applying the operator id^ to ( |6.18| ), gives 

+ 2)d^dACp + ld^'d^da<t>, 



A' 



and, finally, subtracting equation (6.2C) from equation (6.21) yields 

X = ^d^dA<l^. 

From equations ( |S.19 ), ( 6.21 ) and ( 6.22| ), we get the following relation: 

(2 + i){s - 1 - -i')dA<P = l{2 + j')dad^'<l> - l7TAd^'d^db<l>, 



and we also have a field equation which is obtained from equation (6.18) by contraction with the spinor z"" 

= z'^'d^dbcl). 



(6.18) 
(6.19) 
(6.20) 

(6.21) 
(6.22) 
(6.23) 
(6.24) 



For the prese nt pu rposes, equations ( 6.23| ) and ( |6.24 ) are the key equations. Note that equation (6.17) is a consequence 
of equation ( |6.23| ), since the operator (2 + 7') is invertible. To analyze these equations most easily, we henc eforth 
restrict to the case that the quantity s — p is positive. The quantity dA(t> and the right hand side of equation ( 6.23| ) 
are each sums of terms homogeneous of degrees to p — 1 in the variable za'- So with s — p positive, the operator 
(s — 7' — 1) has a well defined inverse acting on such quantities. So we may rewrite equation (3.23) as follows: 



= lis - 1 - 7')-i(a„9^> - (2 + 7')-V^a^'a^afc^ 



(6.25) 



We need to check the com patib ility of equations ( |6.24 ) and ( |6.2| ). First we s how that the Oa derivative of the 
righthand side of equation ( |6.24| ) vanishes modulo the equations ( |6.24 ) and ( 6.25 ): 



dAz'^'d^dbci) = z^'d^dbdAC^ = iz'^'d^dbis - 1 - i)-\dad^'<i> - (2 + 7')- V^^^'a^a,^) 
= i{s - + 7')"' ((1 + i')z'''d^d,dad^'p - z^^'d^dbTTAd^'d^'d^ci)) 

= lis - 7')-i(l + 7')"'(-(l + 7')9f - z^'dad'^'d^dtq^ + ^^e^'^'^f 9b9^5,<, 



= (V2)(s - 7')"'(i + 7')"'(-(i + i')id.)Aacj) - 2z^ da d^dhcj) - 2r 



.J^A'B' 



dad^db 



^it/2)ij'^s)-\id.)AD^ + 2d^ 



B A'B' 



dadb 



(V2)(7' - s)-\id^)An4> + d^eABm = 0. 



Finally we need to show that applying the operator to the righthand side of equation (3.25) gives zero. So it is 
sufficient to show that the quantity i2 + "f')d^dad^ (p + ir^dAd^ d^dbcf) vanishes, modulo equations (6.24) and (3.25). 
We see this as follows: 
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(2 + j')d^dad^'cb + 7r^dAd'''d^db(f> = (3 + y)d^dad^'cl) + d'^'d^dbTl^dACI) 

= id''dA'{s - 1 - 7')"' ((2 + i)dAC'd^'^ - TTAd'^'d^dbcl)) + id'^'d^dbis - 1 - -fT'^^dad^'<l> 



Thus equations (6.24 ) an d ( |6.25|) ar e integrable. If we now write — X]fe=o "^^fc' "^ith 70^ 
get from equations (3.24) and ( |6.25 ) give the foUowing equations: 



= z^'9f4</'fe, 



dA4'. 



fe+i 



i{s-p + k)-' {dad^'^k -{I+P- ky^Ad'^'d^dbcbk). 



k(j)k, for < k < p, we 

(6.26) 
(6.27) 



Equation ( 6.26| ) is vahd for < k < p. Equation ( ^.27 ) is vahd for < A; < p. Note that equation (6.27) entails a 
recursive formula for the quantities <^fe: 

0fc+i = z(fc + l)-\s~p + k)-\z^dad^'c^k - (1 + P - k)-^z^TTAd'''d^db<l>k), for 0<k<p. (6.28) 

Equation ( 6.28|) show s exp licitly that the entire field (f) is uniquely determined by the field 0o- The integrability of 
equations ( 6.24 ) and ( 6.25| ) shows that the system of field equations ( 6.1§| ) for the field (/)(x'', tt"^, z") is completely 
equivalent to the single equation d^db<j)o = 0, for the field 4'q{x°- , n^, za'), in the case s > p. 
Summarizing, we have outlined a proof of the 

Theorem VI. 3 For any pair of positive integers {p, t), there is an isomorphism between the twistor sheaf cohomology 
group H^{U, 0{p, 2p + t — 2)) and the space of holomorphic solutions of the field equation z^ d^db4> = 0, where the 
field (j) is homogeneous of degrees (p + t^p) in the variables (tt^, za')- 

In particular, we have the 

Corollary VI. 4 The space of holomorphic solutions of the equations dA{A'4'B^"c') ~ ^ f'^^ spinor fields with m + 1 
unprimed and m primed indices is isomorphic to the twistor sheaf cohomology group H^{U, 0{m + 1, to)). 

We would like to make several remarks at this point. Firstly, although we have shown the existence of the twistor 
correspondence for several different kinds of fields (i.e., with different index structures when considered as spinor fields 
on space-time) it is only the fields with homegeneity (p + l,p) that can consistently propagate on a curved manifold. 
All other fields suffer from the existence of consistency conditions. The twistor treatment in this work is to some 
extent new in that we use an afhne t wisto r space. This allows to incorporate all homogeneities into one formula (see 
e.g., equation (3.4) in comparison to (S.6)). 



C. The group representation 

We observe that there are natural operators acting on the twistor cohomology groups H^{U,0{p,q)). Indeed 
consider the operators P"^ and Qa/s and E"^, which act on a representative function F{z,(), homogeneous of degrees 
(p, q ~ p) in the twistor variables (z, C) and obeying the differential equation <^ • dzF — 0, as follows: 

P(F) = (z A C)F, (6.29) 

Q(F) = (9, A dc)F, (6.30) 

E{F) = (z ® 9^ + C ® 9c + S)F. (6.31) 



In equation ( |6.31 ), the operator S is the Kronecker delta tensor acting on the representative F by multiplication. 



Note that each of these operators commutes with the operator ( ■ dz- The operator P gives a map from H^{U, 0{p, q j) 
to H^{U, 0{p + l,q + 2)), the operator Q maps H^{U, 0{p, qj) to H^{U, 0{p — l,q — 2)) and the operator E maps 
H^{U,0{p,q)) to itself. These operators generate a Lie algebra under commutation. Indeed, by direct calculation, 
we have the following commutation relations: 

[P"^, Q^s] = -^[>f5', [E], P''"] = -24"P^1^ (6.32) 
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A dimension count gives dimension twenty-eight for this algebra, six for each of the operators P and Q and sixteen 
for the operator E. The operator E generates the complex general linear algebra G'L(4,C). The algebra G'i(4,C) in 
turn is isomorphic to the conformal orthogonal algebra C0(6,C) (the orthogonal algebra together with a dilation). 
Adding in the operators P and Q to this algebra gives the complete algebra of 0(8, C), regarded as the conformal 
algebra associated to 0(6, C), with P forming the translations, Q the generator of special conformal transformations, 
the tracefrec part of E generating rotations and the trace of E giving the dilation. If we introduce the standard 
pseudo-hermitian form on twistor space of signature (2,2), then this algebra has the natural real form 0(4,4), with 
the operator iE self-conjugate and Q the pseudo-hermitian conjugate of P. So we have shown that the direct sum 
over p of all the cohomology groups {U, 0{p,2p + t — 2), gives, for each fixed t, a complex representation of the Lie 
algebra of the group 0(4, 4). It remains an open question whether or not this representation is unitarizable. Indeed 
the "natural" inner product, derived from the action of section one above is not positive definite in the case of spin 
greater than one half. So it would seem that the representation is "naturally" defined on a space with a "natural" inner 
product, but not a Hilbert space. If one took such representations seriously, it would apparently require enlarging the 
framework of quantum mechanics to accomodate " negative probabilities" . 

Finally we note that although this algebra is most easily derived in the twistor picture, one can easily translate 
into the spacetime picture, using the techniques of this section. In the spacetime picture the operator E acts on 
the fields as the (complex) conformal algebra of spacetime. The operators P and Q and E act as follows on a field 
(f){x"','KA' jT^^) obeying the equation M'(j) = and homogeneous of degree s in the spinor ir^: 

P'"^'^ = (^Ti)^''"'''^''''"'^^''^' ^^"^^^ 

= -^^^(5,)c5f XgQ^]C'0, (6.34) 

= -iX^'«X|'5d<^ + X'''"5c'07rn'd^'<p - X^S^{d^)D{7r^cj>), (6.35) 
where we have used the following definitions: 

Pg = {s + 1)6% - iX^'''db, QcB' = {s + l)<5„s' - iX^db, (6.36) 

X^''- = {ix^'^,5%), X^ ^{S^,-zx^^'), 5%^ (6^,0), <5«B' = (0,4). (6.37) 



VII. CONCLUSION 

We have presented the properties of a class of linear equations for fields with half integer spin which generalize 

the Weyl equation for a neutrino. We have shown how the equations arise as Euler-Lagrange equations for a variational 
principle. The equations are of hyperbolic type in the sense that the Cauchy problem is well posed and that there 
exists the notion of a domain of influence. The characteristics of the system are multiply sheeted. The fields propagate 
freely on any curved background, i.e., there are no constraints on spatial hypersurfaces to be satisfied by the Cauchy 
data. The solutions lie in certain Gevrey classes provided that the Cauchy data and the metric of the underlying 
manifold do so. We find a strong relationship between the spin of the fields and the smoothness of the metric, ranging 
from only C*^ in the neutrino case up to analyticity in the limit m oo. We analyzed the characteristic initial value 
problem using the formal method of exact set and showed that it is well posed in the curved background case as well 
as when the system is coupled to gravity via the Einstein equation. It is interesting to note, that this is a system of 
partial differential equations that is not symmetrically hyperbolic (unless m = 0) but still allows the description via 
an exact set. All other examples of exact sets so far have been systems of equations which were also symmetrically 
hyperbolic. This implies that those two characterizations are not mutually included one in the other. We have given 
the general solution of the equations in Minkowski space by first solving the equation □'"+^0 = using Fourier 
methods and then deriving the Fourier representation for positive frequency fields. Finally, we presented a twistor 
correspondence between the the space of holomorphic solutions and sheaf cohomology groups on projective twistor 
space. 

The solution space of the equations in fiat space is a representation space of the Poincarc group. In contrast to the 
case of the massless free fields, however, this representation is reducible unless m = 0. This can be easily seen from 
the fact that the entire solution space for spin 2m — 1 is mapped injectively into the solution space for spin 2m -|- 1 
by the operator L. The solution space is also a representation space for the conformal group. It is not yet known 
whether this representation is irreducible. 

It would be interesting to find similar classes of consistent higher spin equations for integer spin generalizing the 
Maxwell equations. So far attempts have been unsuccesful. 
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APPENDIX A: GEVREY CLASSES OF FUNCTIONS 

Essential in the proof of existence and uniqueness of solutions of non-strictly hyperbolic systems of partial differential 
equations is the notion of Gevrey classes of functions. These are sets of C°°-functions, labelled by a real number a > 1 
which in some way interpolate between analytic functions (a = 1) and functions which are only C'' (conventionally 
made to correspond to the case a — oo, see below). 

Deflnition A.l Let S be an open set m K', p > 1 and a > 1. Then 7p"''(5) is the set of functions f : S —> C such 
that 

sup ^1 \\D^f,S\\t (Al) 

where a is a multi-index a = (cri, . . . , ai), |cr| = (Ti + . . . + ctj and ||/, 511^ is the usual U -norm of f . 
Similar classes are defined to characterise the behaviour of the functions with time. 

Definition A.2 Let S := [0,r] x 5 &e a stri-p in M'+V. -p > \, n > \ and a > \. Then 7p'^"^(E) is the set of all 

functions f : —i- <C such that 

where, again, a and (3 are multi indices with (Tq = (i.e., a refers only to "spatial" derivatives) and < x'^ < T. St 
is the slice x'^ = t, as usual. 

We extend the definition with a = oo by the rule that (1 + = 1 for \a\ — and (1 + \(t\)" — otherwise. Then 
lp°°\S) is equal to the function space Lp{S) and 7p'^°°''(E) is equal to the set of all functions for which \\D^ f, St\\ is 
a bounded function of t for all /? with |/3| < n. 

For p = oo, I = 1 and for real valued functions this is the classical case of Gevrey |^. In that case j^(^\s) is an 
algebra which is closed under the composition of its elements. A similar property holds in the general case (see |^). 

Here are some basic properties of the Gevrey classes: they grow with a; if ai < a2 then jp^^"'\i:) C Jp^^"'\e). 
Also 7p'^"'*(E) C 7™'^"''(E) if m > n. If a = 1 the classes consist of functions which are analytic in x^, . . . , x', but for 
a 7^ 1 one can show that there exists a partition of unity into elements of the classes with arbitrarily small support; 
functions with compact support are not necessarily zero. 

The essential qualitative distinction within the Gevrey classes seems to be between the case a = 1 and a > 1, 
the latter case permitting domains of influence and thus allowing the study of wave propagation. This indicates the 
hyperbolic character of the equations under consideration. The second distinction is between the cases of finite a and 
a = oo. Infinite a permits the existence of only a finite number of derivatives and thus the appearance of shocks is 
possible indicating the strictly hyperbolic case. 
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